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Optical Absorption Lines of Hydrated Manganous Salts 
at Low Temperatures} 


By R. PappaLarpo 
H. H. Wills Physics Laboratory, Bristol 


[Received July 18, 1957] 


ABSTRACT 

The absorption spectra of four manganous hydrated salts (sulphate, 
chloride, fluosilicate and ammonium sulphate) have been investigated at 
78°K and 20°K in the 15 000 cm! to 30 000 cm-! range. 

Several groups of lines are found, some of them showing subsplitting of 
the order of 10 em~1. 

The results are satisfactorily interpreted by means of the ligand field 
theory ; a qualitative discussion of refinements of the cubic field scheme 
in relation to the observed spectra is also given. 


§ 1. INTRODUCTION 


THE absorption spectra of complex ions of elements of the iron group 
have been thoroughly investigated in recent years for the case of ions in 
solution (in particular, by Jorgensen 1956). Though the low temperature 
optical absorption has been studied in crystals of many salts of rare-earth 
elements, for the iron group (with the exception of Cr*+) there is only the 
preliminary work of Gielessen (1935) and the recent work of Holmes and 
McClure (1957). The existence of a rather crude but successful theoretical 
scheme for interpreting optical absorption (Tanabe and Sugano 1954, 
Orgel 1952), and the achievements and limitations (Abragam and Pryce 
1951 a, b) of the paramagnetic resonance methods, have both helped to 
stimulate considerable interest in the subject: it would be worth while 
making a systematic study of the optical absorption for the iron group 
of elements. 

The theoretical scheme of interpretation—common to the cases of 
solutions and crystals—is based on the ligand field theory (‘Tanabe and 
Sugano 1954, Orgel 1952). The complex ion is considered as a system of 
four or six ligands bound to the centrai metallic ion by polarization forces ; 
this produces an electric field of determinate symmetry, with which the 
3d electrons of the iron group element interact. Hence the spectroscopic 
terms of the free ion are altered, and the optical absorption is due to 
electron transitions between 3d orbitals, giving the crystal its charac- 
teristic coloration. These transitions, forbidden in the free-ion state, 
have a low probability in the complex ion. 


+ Communicated by Professor M. H. L. Pryce, F.R.S. 
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For the theoretical approach to the problem it is necessary to calculate 
the new set of levels of the ion; for this purpose the ligand field model 
uses pure ionic functions which must possess the symmetry properties of 
the ligand field. The group theory treatment by Bethe (1929) gives 
readily the qualitative trend of the levels, in the cases of weak, inter- 
mediate, or strong fields with respect to the spin-orbit and spin-spin 
interactions. This calculation applies to crystals in general. Quanti- 
tative data for the resulting levels of the iron-group elements have been 
calculated by Tanabe and Sugano (1954) for the case of strong field, 
cubic symmetry (six octahedrally coordinated ligands), and neglecting 
spin-orbit interaction and non-cubic components of the ligand field. 
Orgel (1952) gave the analogous treatment for a weak field. 

This formalism has been widely used to interpret the absorption bands 
of solutions; the agreement between theory and experiment is quite 
satisfactory. A more exacting test of the ligand field theory would come 
from the study of crystals, since by working at low temperatures it is 
possible to resolve details of the absorption spectra, masked at room 
temperature by the broadening due to thermal vibrations. The study of 
the optical behaviour of crystals at low temperatures would give infor- 
mation on the departure of the ligand field from cubic symmetry and on 
the roles played by spin-orbit coupling, configuration interaction, and 
lattice vibrations. Information on the actual symmetry of the field 
acting on the 3d electrons in complex ions, and more accurate knowledge 
of their ground and excited states, would be useful for the interpretation 
of magnetochemical and paramagnetic resonance data (Abragam and 
Pryce 1951 a, b). 

An experimental study of the absorption spectra of hydrated crystals 
containing iron-group elements has therefore been initiated. The present 
paper reports results on various manganous salts. 


§ 2. EXPERIMENTAL 
2.1. The Cryostat 

The basic requirements of the cryostat were to cool the crystal and 
keep its temperature constant for a reasonable time, to avoid formation 
of frost along the path of light, and avoid the passage of light through 
the refrigerating liquid. All these requirements were met by the cryostat 
of fig. 1; the main disadvantage was that only one crystal could be 
investigated for each filling of hydrogen or nitrogen. 

The cryostat was designed by Dr. H. D, Keith, and in its main features 
was similar to the one used by Johnson and Studer (1951). The trans- 
mission cell and the container for the refrigerating liquid were made from 
a single copper block so that if the crystalline sample was in good thermal 
contact with the body of the transmission cell it would reach the tem- 
perature of the liquid bath. To maintain constant temperature in the 
crystal for a reasonable time, the transmission cell and inner container 
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Fig. 1 


w> SZ 


1O5¢m : 
The cryostat. 


A, inlet for refrigerating liquid ; B, hydrogen gas outlet ; ©, exchange gas inlet ; 
D, greased joint; E, liquid nitrogen inlet; F, corrugations in glass 
jacket ; G, copper tubing for exchange gas ; H, main pumping tap for 
interspace vacuum; I, glass portion of liquid hydrogen container ; 
L, Kovar seal ; M, main copper block ; N, transmission cell ; O, waxed 
windows ; P, liquid nitrogen. 
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were surrounded by a vacuum space and enclosed in a second dewar filled 
with liquid nitrogen. The nitrogen jacket was designed so as not to 
intercept the beam of light passing through the sample ; and this helped 
to maintain the outermost windows of the cryostat at room temperature 
free of condensation and frost. 


2.2. The Transmission Cell (fig. 2) 

This proved to be the most troublesome part of the system. It was 
originally designed for work with liquid helium; since at this tem- 
perature good thermal contact is essential, the cavity contaiming the 
sample was filled with helium exchange gas at low pressure. To main- 
tain the cavity gas-tight at low temperatures the following method was 
used. Transparent windows were mounted in thin copper diaphragms 
with Araldite and baked for an hour at 190°c (Roberts 1954); the 


Fig. 2 


32cm 


The transmission cell. 
A, Hoyt metal soldering ; B, thin copper diaphragm ; ©, transparent window ; 
D, Araldite ; E, small hole for gas exchange circulation ; F, crystalline 
sample ; G, copper ring; H, light spring ; I, small screw to lock the 


nai a ring L; M, exchange gas inlet; N, tapered hole for the light 


diaphragms were then soldered to the main copper block with Hoyt 
metal No. 314. One of the windows was left assembled permanently 
and the other could be removed and replaced when changing the sam i 
The heat required for soldering was liable to spoil the crystal; on ‘the 
other hand leaks would sometimes occur in the soldering during coolin 
spoiling the interspace vacuum and causing rapid evaporation of the 
refrigerating liquid. 
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For the work with liquid hydrogen, the thermal contact obtained by 
clamping the crystal was quite satisfactory. The crystal was placed in 
the transmission cell in front of a tapered aperture, and pressed against 
the main copper block by a copper ring and a weak spring. It was locked 


in position by another copper ring, which was also i ack ; 
with the body of the cell. Page cou 


2.3. The Spectrograph 


For observing the narrow absorption lines, the more sensitive photo- 
graphic method was preferred to spectrophotometric apparatus employing 
photocells and a monochronomator. A Hilger E.1, large quartz spectro- 
graph, was used to investigate the range from 33004 to 64004 
(30 300 cm™! to 15 600 cm). Typical values of the dispersion in various 
spectral regions were as follows : 


7A/mm (60cm-Y/mm) at 34004 (29400 cm~}), 
15A/mm (70cm-4/mm) at 42004 (23 800 cm~}), 
(100 cm~1/mm) at 50004 (20000 cm~}), 

( (16700 cm-*). 


25 A/mm 


35 A/mm (100 cm~4/mm) at 60004 


The light source was a 50 watt lamp with a vertical filament. The light 
was first focused on the crystal sample and then concentrated on the slit 
of the spectrograph, using two quartz lenses. The spectrograms were 
taken on Ilford Rapid Process Panchromatic plates, and developed in a 
high contrast developer. On each plate some 12—i4 spectrograms were 
taken, with exposure times increasing from 1 second to 10 minutes. The 
line spectrum of a mercury vapour lamp was recorded on each plate for 
calibration purposes. 


2.4. Wavelength and Intensity Determinations 


For wavelength measurements a precise plate calibration was derived 
using an iron arc spectrum. Reference points on this curve were given 
by the mercury lines recorded on every plate. All wavelength measure- 
ments were made by scanning the plates with a Hilger microdensitometer ; 
photographic enlargement was used as a check on the photometric tracing 
in any case of doubt. 

Intensity determinations were restricted to absorption lines, since an 
accurate determination for bands is not possible with the photographic 
method. To evaluate the intensity of a line a photometric trace of the 
plate was taken, choosing a suitable exposure time. From this the 
general behaviour of the absorption near the line was deduced. To 
correlate quantitatively the value of the absorption peak of the line to its 
extrapolated background, it was necessary to know the sensitivity curve 
of the plate in that particular region, and for the amount of light passing 


through the given sample. 
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For this purpose galvanometer readings were taken, measuring the 
blackening of the plate corresponding to various exposure times for a 
given wavelength ; for good accuracy the setting was made at 1 or 2mm 
from the line. A set of curves (galvanometer deflection against exposure 
time) was thus drawn ; on one or more of these curves, points were located 
corresponding to blackening at the maximum of the line and to the 
blackening of the extrapolated background in the absence of lines. Since 
far from cut-off the sensitivity of the plate did not vary appreciably over 
1 or 2 mm, it was reasonable to assume that the deduced ratio of exposure 
times for the extrapolated background and the peak of the absorption 
line, was very nearly the value of the ratio at the actual site of the line. 
This ratio of exposure times was then correlated with the corresponding 
intensity ratio. To minimize errors, only the linear part of the cali- 
bration curves, far from saturation, was taken into consideration. 

The values of the absorption deduced in this way did not require 
correction for reflection losses, were not sensitive to slight alterations in 
the geometry of the system, and did not require constant developing and 
fixing conditions from plate to plate. Naturally the procedure could not 
be applied to very strong absorption lines. 


2.5. The Crystalline Samples 


The crystals were grown from solutions of analytically pure chemicals. 
The samples were about 44mm across and had thicknesses varying 
from 1mm to 7mm. They were lightly polished to reduce reflection 
losses and interference effects. 

Samples of MnsSiF,.6H,O and Mn(NH,).(SO,),.6H,O were kindly 
supplied by Dr. L. C. Jackson. 


2.6. Operation 

The central part of the cryostat (liquid container—transmission cell) 
was demountable. After inserting the sample, this removable part was 
replaced and the cryostat sealed on top by a greased joint. The inter- 
spaces A and B and the transmission cell were evacuated. (Loss of 
transparency during this operation and the subsequent cooling sometimes 
spoilt the crystal samples.) The interspaces were then filled with dry 
oxygen gas at 0-5 atmospheric pressure, to provide an agent for heat 
exchange in the preliminary cooling of the system. The inner liquid 
container was also evacuated and flushed several times with hydrogen gas. 

The outside jacket was slowly filled with liquid nitrogen. When the 
transmission cell reached the temperature of liquid nitrogen—usually 
after about 1} hours—the interspaces were again exhausted, and liquid 
hydrogen (or nitrogen) was syphoned in. 

Usually the 250-300 cm of liquid hydrogen lasted for 5-6 hours ; in 


this time several spectrograms could be taken, once the temperature of 
the system had stabilized. 
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§ 3. EXPERIMENTAL ResuLts 


Absorption spectra of each sample were taken at room temperature, at 
78°K and 20°K. Only data for the measurements at 20°K are tabulated 
below ; an outline of the results at the other temperatures is also given, as 
this is useful for the tentative interpretation that follows. . 


Fig. 8 


Gaivanometer 


deflections 


(arbitrary units) 


Microdensitometer curve for MnSO,.4H,O at 20°K (25 000 cm-! region). 


The formula (3.1) of Tanabe and Sugano (1954) was used for the deter- 
mination of oscillator strengths. The tabulated extinction coefficient for 
the maximum of a line is the coefficient « of the usual relation 

Leal, Lise: 
The densities, molecular volumes and crystal structures of the salts were 
taken from Mellor, T'reatise on Inorganic and Theoretical Chemistry. As 
regards the accuracy of the tabulated data, there may be occasional dis- 
erepancies of 80-100 cm~! in the absolute wave-numbers : differences in 
wave-number of more than 100 cm! should be fairly accurate. The 
aecuracy of the values of the absorption coefficient is not very good, and 
only the first figure is significant; this is reasonable considering the 
limitations of the photographic method. Values of relative intensity 
are only given for lines which have similar intensities and are not widely 
separated. 
3.1. Manganous Sulphate MnSO, . 4H,0 

The erystal is monoclinic, with axial ratios a:b: c¢:: 0-8643: 1: 0-587, and 

B=90° 53’. The sample investigated at 20°K (fig. 5, PI. 56) had a thickness 
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of 695mm. In the spectral region investigated (3300 A—6400 A) there 
were several groups of lines (see table 1). 


Table 1. Lines of MnSO,.4H,O at 20°K 


: Sub-splitting 
| Group | Wave- Half- ; Extinct. ; : A “ 
of enor width | Relative coeff. Oscillation | = a] Remarks 
Tinea” Panisdee in intensity max strength = Bo 
7 em} (em—?) & 
24838 30 0-5 0-3 Diffuse 
A 24882 30 0-6 0-3 alms band at 
RT 
24950* 30 il 0-4 Ipc lOm es Present at 
RT 
25025 weak 
25100 65 0-45 0-1 0-8 x 10-8 
25176 37 1- 0:3 Hell Se 2 ~15 | rather 
broad 
B 25221 0-65 2 ~25 | rather 
ee 68 0:8 0-2 isGyse UO possibly broad 
2 0:35 
25826 | ~30 ~3-01 B54 I) weak 
Cc 25880 2 ~30 | very weak 
but sharp 
25960 0-03 weak 
D 27292 27 0-14 0-4 10-8 
27352 28 0-12 0-3 x 10-8 
see 54 0-12 0-66 x 10-8 
27640 
211ol 0-4 
27778 
27855 
27910 
27980 0:3 
28050 
E 28090 
28177 
28248 
28280 
28296 
28369 
See 48 0-02 | 13x 10-° 
28498 
28670 
F 29700 


very diffuse 

weak 

_ Se eo 
Oscillator strength 1-4 x 10-8 for the same line at room temperature. 


A. A rather intense triplet centred on 24 900 em-! (fi 

‘ent ; 6, BIT 56 3 fi 
B. A complex pattern of lines at 25 200 cm-! 1 ge 6, Pl. 56; fie 3) 
C. A faint quadruplet at 25 800 cm! (fig. 6, Pl. 56: fig. 8) 
ie a intense doublet at 27 300 cm—? (fig. 7, Pl. 56; fig. 9) 


. A broad band of complex structure centred on 27 800 cm—! (fi 
Xs 7 
fig. 9). The peaks of this band are shown in table 1. eT aloe 
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) 3.2. Temperature Effects 

At room temperature (fig. 3, Pl. 55) group A consists of a line at 
24 934 cm™ and a hazy band at 24 820 cm-1: group B is a diffuse band 
from which two lines can be distinguished (at 25 120 cm-1 and 25 260 cm~?). 
At nitrogen temperature (fig. 4, Pl. 55) both groups A and B are already 
assuming their final form—subsequent cooling to 20°x (figs. 5-7, Pl. 56) 


Fig. 9 


Galvanometer 


detle<at tons 


Carbitrary units ) 


Microdensitometer curve for MnSO, .4H,O at 20°K (27 000 cm region). 


sharpens group A and reveals the subsplitting in the lines of group B. 
Group C only appears at 20°K. Groups D and E are absent at room 
temperature, and there is only a very faint, fairly broad band. Doublet D 
appears at liquid nitrogen temperature, together with a few weak lines 
of group E. At 20°K the doublet D is sharpened, and the complex pattern 


of E becomes very clear. 


3.3. Manganous Chloride MnCl, . 4H,O 


The crystal is monoclinic, with axial ratios a@:6:c¢:: 1:15: 1: 0-64, and 
B=99° 25’. The thickness of the sample investigated was 1-37 mm. 
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Table 2. Lines of MnCl, . 4H,0 


2 
oo 
DQ 
= 
=e 
ch. 
5 
itle} 


Wave- | Half- Extinct. 


Group number | width | Relative coeff, Oscillator 2 os Remarks | 
aS a in intensity max strength = g 
lines 20°K em-! (em?) s 

24318 26 0-7 0-3 10-8 | 2 ~ 9cm™ | sub- 

: splitting 
very 
sharp 

A 24462 45 0-6 0-3 1-5x10-® | 2 ~ll cms 
24500 42 0-35 : 
24594 42 1 0:3 1-7x10-§ | 2 ~ 6em™ only line 
possibly at room 
temp. 
24690 0:3 ae he 
B 24802 45 0-7 : 

24919 60 1 ee lO ae 2 possibly 

25300 

27137 1) 1 0:8 1-8 x 10-8 sharp 

C 27174 26 0:97 equi- 

27207 28 0-96 J ba 

ines 

27300 

27426 

D 27525 fa 1 0:5 4-5 x 10-8 possibly 

27570 1 

27770 

29070 Te 1 1 OSs! 3 possibly 

E 29155 95 1 3 possibly diffuse 

29330 


At 20°K the absorption spectrum is composed of (see table 2) 


A. A main triplet with various subsplittings, at 24 500 em—! (fig. 12, Pl. 57; 
fig. 14). 

B. Faint lines at 24 800 cm—! (fig. 14). 

C. A strong triplet of sharp lines at 27 150 em! (fig. 13, Pl. 58; fig. 15). 

D. A band extending up to 28000 cm-!, showing some structure. The 
peaks of this structure are shown in table 2 (fig. 13, Pl. 58; fig. 15). 


E. Two broad lines at 29100 cm-!, showing some structure which is not 
very definite even at 20°x (fig. 13, Pl. 58; fig. 16). 


3.4. Temperature Effects 


At room temperature group A is very faint, but not diffuse, and the. 
line at 24 600 cm~ is the strongest of the group. The whole pattern 
becomes stronger at liquid nitrogen temperature (fig. 10, Pl. 57) ; finally at 
20°K the subsplitting is enhanced, and the lines of group B too become: 
more distinct (figs. 11-13, Pls. 57, 58). The other groups are absent at room 


temperature. They first appear at liquid nitrogen temperature, though. 
rather blurred, becoming much sharper at 20°K. 


ee a ee 


w=" - | Ss ee SC 
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Galvanometer 
deflections 
(arbitrary units) 


20 


24300 = 24600 24-900 25200 em! 


Fig. 15 


20 


Galvanometer 
deflections 
@rbitrary units) 


15 


27-000 27-400 


+ _—_____—_ 


Microdensitometer curve for MnCl, . 4H,0 (~27 000 cm-'), 
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Galvanometer 


deflections 


1S Carbltrary units ) 


29:6 x10? 


Microdensitometer curve for MnCl, .4H,O at 20°K (29 000 em-? region). 


Fig. 18 


Galvanometer 
deflections 
1S 


Carb Units) 


— = aie att ol 
25 25°3 25°6 259 x10? ecm"! 


Microdensitometer curve for MnSiF, . 6H,O at 20°x. 
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3.5. Manganous Fluosilicate MnSiF, . 6H,0 


The crystals belong to the trigonal system, and have the axial ratio. 
@:c¢::1:0-504, and «=112° 30’. The thickness of the sample used was 


3°87 mm. 
Table 3. MnSiF, . 6H,O 


rel Wave- | Half- Extinct. Subsplitting 
roup : d ‘ : = 
of number width Relative coeff. Oscillator 3 ey 
fea : a c in intensity max strength = 2 E Remarks | 
20°K em (em-!) 3 
5 
24890 _Wealk 
sharp 
doublet 
A 24946 
25002 2 & 
25075 65 0-25 Seales | Present, | 
25170 rather |; 
blurred 
at room 
B 25330 85 0-45 7x10-8 | 9 ~14 | | temper- 
25420 | 65 0-5 6x 10-8 SE a 
Cc 26413 2 
28900 ? 


At 20°K the spectrum consists of (see table 3 and fig. 17, Pl. 58) 


A. One intense line, sided by a faint doublet at 24 980 cm-! and a broad 
line at about 25 170 cm-}?. 

B. A fairly intense triplet, with two of the lines close together. 

C. Very weak and narrow bands at about 26410 cm and 28 900 cm-! 
though the evidence from the photometric tracing is not very conclusive. 


3.6. Temperature Effect 

As for the other salts, the cooling clarifies the relatively strong lines. 
already existing at room temperature, and subsplitting becomes apparent 
at 20°K. 

3.7. Manganous Ammonium Sulphate Mn(NH,).(SO,4). . 6H,0 

The crystals have axial ratios a: b: ¢ : : 0-736: 1: 0-497, and B=107° 2’. 
The sample used had a thickness of 1:98 mm. Only one line is detectable 
at 78°k, its wave-number being 25 220 cm~'. Cooling to 20°K produces. 
no appreciable change in position or intensity. At 20°K the line is pre- 
ceded by a narrow band centred on 25 050 cm™' and by a very weak 
sharp doublet at 24 750 cm~? (fig. 19). Samples of (Mn, Zn)SO, with 
dilution ratios 1:1 and 1:5 were prepared to study the effect of 
diluting the paramagnetic ion. Unfortunately, both samples were spoilt 


during the cooling. 
§ 4. INTERPRETATION OF THE EXPERIMENTAL RESULTS 
4.1, Oubic Field Model 


We shall now compare the theoretical predictions with the experimental! 
data for the Mn2+ ion. Firstly it must be remembered that the caleu- 
lations by Tanabe and Sugano were made with the hypothesis of a cubic 
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ligand field, and that distortions from cubic symmetry and spin-orbit 
effects were ignored. Consequently, from this model we can expect to 
have an explanation of only the main features of the absorption. 


Table 4. Mn(NH,),(SO,)2 . 64,0 


ee 


rel Wave- | Half ; Extinct 
ee P | number width Relative coeff. Oscillator Sub-splitting | Remarks 
2 at in intensity max strength 
ules 20°K emp? (Gr) 
24728 ee 
24795 | ~15 ell Ome sharp 
J lines 


The predicted levels for 3d systems were plotted by Tanabe and Sugano 
as a function of two adjustable parameters, 4 the ligand field strength, 
and B a combination of Slater-Condon parameters. For the Mn?* ion 
the authors assume values of 1230 cm~! for 4, and 860 cm~! for B, as for 
the free ion. However, the agreement with the absorption spectra of 
solutions is rather poor ; Orgel (1955) suggested a value of 780 cm~! for 
B, and then there is better agreement between the theoretical levels and 
the observed bands in solutions containing [Mn(H,0),|?*. It should be 
noted that this decrease in the value of B from the free ion to the complex 
ion state seems to be a common feature in other elements of the iron group, 
especially for trivalent ions ; this is interpreted either as a spreading of 
the electronic cloud of the ion, or as evidence of partial covalent bonding. 
Assuming Orgel’s value for B, we have the following sequence of quartet 
states beyond the °A,(de?, dy?) ground state. 

47’, (de*, dy?) at 15 200 em-1} 

4T, (de*, dy) at 19 500 em—1y4 

4A, *Ea (de? dy?) at 25 300 cm—! 

“T, (de®, dy?) at 26 950 cm~! (extrapolated value) 

4B (de3, dy?) at 30 400 cm—! 

*T, (de®, dy?) at 31 200 em~! (extrapolated value). 
Levels having the same subshell configuration (de>Ndy®) possess the 
same cubic field dependence, so that transitions between such states are 
expected to produce sharp lines instead of the more common bands. 
Transitions *A, (de®, dy?) > 4T,, *T,(de4, dy) are expected to give only 
weak bands. Before applying these predictions to the crystals investi- 
gated, we must consider the evidence that justifies the assumption that 
the Mn** ion is surrounded by six. water molecules, located approxi- 
mately at the vertices of an octahedron. From x-ray analysis it is known 
that this is certainly the case for MnSiF,.6H,O. This has a distorted 
CsCl structure of octahedral radicals of Mn(H,O), and SiF,, of similar 
ee ae ee 


+ In aqueous solution Jorgensen (1956) finds these levels at 18800 and 
23 000 cm-! respectively. 
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shape and size, and lying parallel. There are no x-ray data for the other 
salts, but it seems likely from crystallographic evidence on all the salts of 
these transition elements that we can safely assume the octahedral model 
as a first approximation since the characteristics of the absorption are 
broadly similar for all the salts investigated. 


4.2. Intersystem Combinations 
In the free ion state for Mn?+ there are no sextet states up to 50 000 cm~! 
apart from the ground state. Hence it is reasonable to consider the 
observed lines in the complex ion as due to intersystem combinations of 
SA, (ground state) with upper quartet states. Transitions to levels of 


Fig. 19 


Galvanometer 
deflections 


(arbitrary wnits) 


Microdensitometer curve of Mn(NH,),(SO,), .6H,0 at 20°K. 


this kind are forbidden even in complex ions (see Hellwege 1948) rt ne 
spin S is a true quantum number. However, if there is a eee ) 
sextet with quartet states in the ground state, these ee ey 
possible even if their probability is small. Since the pronsy i of the 
transitions is proportional to the square of the ee in i e bs 
functions, assuming an oscillator strength ii of a for a ar oS 
combinations and since the observed value of fis ~10-8 for the et a ce 
manganous ion, the amount of contamination would be of ute a er We 
which is reasonable. intersystem combinations with onc pone : 
improbable because the spin-orbit coupling couples leve Me : bates 
Even if possible, they would only give very weak bands and no 
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3d systems. If we consider the spectra of MnCl, . 4H,0 and MnSO,. 4H,0 
we find that there are lines at ,~25 000 cm~!; ~27 000 cm! and 
~29000 em-1, as expected if transitions occur from °A,. to 4A,, “Ka 5 
4T, and: ‘Hy respectively. Moreover, the wave-numbers of these lines 
are independent of temperature (even in the wide range from room 
temperature down to 20°K), as predicted. These are certainly remarkable 
achievements of the theory. 
4.3. Selection Rules 

Groups of lines appearing in some salts are missing in others; this 
could be explained in terms of selection rules that are dependent on the 
crystalline symmetry. Hellwege has calculated a set of selection rules 
valid for crystals of various symmetries, for electric and magnetic dipole, 
and electric quadrupole transitions. These selection rules cannot easily 
be applied in the present case since neither the symmetry and strength of 
the non-cubic part of the field, the order of magnitude of spin-orbit 
coupling, nor the nature of the transitions involved (electric or magnetic- 
dipole, or electric quadrupole) are known. Moreover Hellwege did: not 
consider the effect of lattice vibrations. 


4.4, Hacited Lattice Vibrations 
In a field of On symmetry electric dipole transitions are not possible 
because the overall parity cannot be conserved in the process. Hence in 
this case the parity requirements can only be satisfied by a combined 
process of electron transition and the excitation of lattice vibrations 
of suitable symmetry. However, this kind of process is not strictly 
necessary in the case of lower crystalline symmetry. 
Some of the observed lines do seem to be coupled with lattice or 
complex ion vibrations, but others do not, although their oscillator 
strength is compatible with that of electric dipole transitions. 


4.5. Spin—Orbit Coupling 

For half-filled shells, as the d>-configuration of Mn?+ and Fe®+ the 
Hund rule gives a spin-orbit coupling which is zero in a first order approxi- 
mation. Second order effects in the free ion give splittings of about 
50-100 cm~*. | We will assume that the spin-orbit coupling is present 
and effective to a similar extent in resolving the spin-degeneracy of 
orbital levels, even if we can only estimate the order of magnitude of its 

ct. 

at 4.6. Comparison of Spectra of Various Salts 

(a) In general for the four substances investigated the number of 
spectral lines increases from the fluosilicate and the Tutton salt to the 
chloride and the sulphate ; this is in order:of lowering of the internal 
symmetry, which is nearly trigonal in the first two cases (from magnetic 
susceptibility measurements) and rhombic in the latter two (as inferred 
from the macroscopic.symmetry properties). 

(6) Consider first the.groups of lines A and B of the various salts. 
From densitometric readings a coupling with vibrations seems, to. be 


i. 
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excluded—at least, in the sense that there is no excitation of a complex 
vibration pattern accompanying the electronic transitions. The influence 
of lattice vibrations seems to be limited to the production of fluctuations 
in the non-cubic part of the field, thus making the lines rather blurred at 
room temperature and at 78°K. Cooling reduces the thermal vibrations 
and thus stabilizes the non-cubic part of the field, the part responsible for 
the separation of the lines; the lines become sharper and small sub- 
splittings of the order 10 cm~! are observed. Now, if our assignment is 
correct—namely, that the lines at 25 000 em! are due to the excitation 
SA, — 4A,, *Ea even in a rhombic field there would be a maximum of three 
lines. Hence it is necessary to postulate that the spin-orbit coupling 
removes the spin degeneracy of the upper levels: the ground state, derived 
from a °S state, for present purposes is to be considered as not split, 
because eventual splittings due to spin-spin coupling are expected to be 
less than 1 em~'in magnitude. On the other hand, for the upper quartets 
the orbital levels are split into sublevels whose characters and numbers 
can be deduced from a group theory treatment (Bethe 1929, Opechowsky 
1940) once the symmetry of the non-cubic part of the field is known. The 
maximum number of levels obtainable from the +A,, 4E level due to spin— 
orbit coupling is six, since all levels remain doubly degenerate by Kramers’ 
theorem. 

We find a maximum of eleven lines in the 25 000 cm~! region of the 
sulphate. Some of these presumably correspond to vibrational excitation. 
Since generally some of the lines found in this region are stronger than 
others, it is tempting to consider that the more intense lines are derived 
from the 4A, level, which for a field of cubic symmetry would be favoured 
in the contamination with the ground state ®A,. 

Concerning the nature of the lines, the oscillator strength of 1078 
observed for most of the lines could equaily well correspond to the theo- 
retical prediction (Tanabe and Sugano 1954) of 10-7 for electric dipole 
transitions, or to 4x 10-® for magnetic dipole transitions, since the pre- 
dictions have an uncertainty of ~10-?. Only a study of optical aniso- 
tropy could resolve the choice. / : oe 

(c) The lines of the quartet found for MnSO, 4H,O at 25 900 cm 
are interesting. They are not a vibration replica of the first strong 
triplet at 24 900 cm~!; they are quite sharp, but too weak to be penned 
with an electric dipole transition (oscillator strength being ~3 X 10-1), 
The lines are probably due to a magnetic dipole or electric quadrupole 
* (a) Now consider the 6A, — ‘F, transition. From the densitometric 
analysis in the 27 000 cm? region it is clear that for the chloride and the 
sulphate the transition is coupled with excitation of vibration modes in the 
In agreement with this, these lines are miore ee by 
changes in temperature than the transitions at 25 000 cm for yan 
lattice vibrations simply produce fluctuations in the field. ; Moreover, 
the lines at 27 000 cm™! are missing in the fluosilicate and the Tutton salt. 


upper state. 
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Now if the field symmetry is higher in these two salts, then correspondingly 
the number of allowed modes of vibrations is less. Hence the transition 
probability for these salts would be severely decreased. 

(e) Finally, the lines observed at 29 000 em! for MnCl, and the 
diffuse band at 29 700 cm-! for MnSO, can be ascribed to the transition 
6A, (de3, dy?) > *H(de?, dy) which is expected at 30 400 em-1}, Since 
the lines of the chloride do not show a vibration structure, in the absence 
of spin-orbit coupling at most two lines are expected, even in a rhombic 
field. The spin-orbit ‘effect increases the number of configurations for 
the E-orbital to a maximum of four. 


§ 5. CONCLUSION 


The ligand field theory applied to Mn?* ions is in general agreement 
with the experimental data. The theory successfully predicts the approxi- 
mate positions of the levels, the character of the lines corresponding to 
transitions to these levels from the ground state, and the fact that the 
wave-numbers of these lines are not affected by changes in temperature. 
However for a detailed investigation of the spectrum, this approach has 
to be refined to include the effects of non-cubic field, spin-orbit inter- 
action and the lattice and complex ion vibrations. The analysis of the 
optical behaviour of complex ions does not appear to be a field for broad 
generalizations. The interaction of complex ions with the ligand and 
radiation fields could only be subjected to a more rigorous analysis by 
studying the optical anisotropy in crystals of high symmetry and the 
optical behaviour in a magnetic field. 
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SUMMARY 


The temperature dependence of the specific heat of simple crystals at 
temperatures well above the Debye temperature is estimated by using, as 
a simple model, the two-atom anharmonic oscillator. It is found that the 
difference between Cy and the Dulong and Petit value is not much less 
than the difference between Cy and Cy. There is rough agreement with 
experiment. 


§ 1. INTRODUCTION 


THE specific heat of solids at constant volume (Cy), derived through 
a well-known thermodynamic formula from the measured value of the 
specific heat at constant pressure (Cp), is normally observed to increase 
steadily above the Dulong and Petit value as the temperature is raised. 
Born and Brody (1921) outlined a quantum theory relating this increase 
to the fact that the forces between atoms do not obey Hooke’s law, and 
a simpler classical treatment was given by Schrédinger (1922). Neither 
treatment makes a quantitative comparison between theory and experi- 
ment. Nor do simplified accounts, using a two-atom model, given in 
books by Kittel (1953), Peierls (1955) and others. 

The following treatment of the two-atom model connects the temper- 
ature coefficient of Cy with Griineisen’s y and the coefficient of thermal 
expansion, and thus allows its magnitude to be estimated. 


§ 2. THEORY 


The standard treatment of the two-atom model considers an atom to 
be executing its thermal motion along a straight line in an asymmetrical 
potential provided by a fixed atom. The total energy F is expressed 
in terms of the momentum p and the displacement x from the equilibrium 


inter-atomic distance (7) by the equation 
q 2 ’ 
B= = +40? + bu*+out+des+.... 
2m 


Here m is the mass of the atom and a, 8, ¢, d, ete. are coefficients deter- 
mined by the shape of the inter-atomic potential. It is important to 


pall EE a ae 
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note that if the equilibrium position is altered by imposing an external 
force on the atom, these coefficients are altered, because the expansion 
of the inter-atomic potential into a power series must then be made about 
a different point of the potential curve. 

The partition function (sum over states) is, in the form given by 
classical statistical mechanics, 


Ge | | exp (—E/kT’) da dp 


where 7’ is the absolute temperature and k is Boltzmann’s constant. It 
is assumed that the effect of the third and subsequent terms of the above 
expansion for # is small, and Z is obtained in the form of a power series 
in kT of which only the first two terms are retained : 


kT 156? se 
om ES (aa) ef 


Here v is the frequency of vibrations of small amplitude and is given by 


l/a\ie 
— = (=) . . . : . . . . (1) 
The heat capacity per atom at constant volume is, to the first order in 7’, 
C2 
Oy=T ar (kT In Z)=k(1-+-2BET) Se eee 
where 
_ 156% = 8e 


This is the usual result. The purpose of the following analysis is to 
relate B to measured quantities. ; 


The mean value of x is (Kittel 1953) 


(v)=Z-1 | exp (—B)kT) de dp=— aD meee ACS 
to a first approximation. It is evident that if |b| is large both the 
thermal expansion and the increase of Cy will be large. This correlation 
noted by Born and Brody, suggests that B should be evaluated by rene 
it to the thermal expansion, and this may be done by introducing also the 
change of vibration frequency with volume. 

A change of volume, in terms of this model, implies that the position 
of equilibrium has been shifted by an impressed force. As previously 
noted, the coefficients a,b, c, etc. are functions of the equilibrium inter- 
atomic distance 7. If the distance from the new position of equilibrium 
is z, and the impressed force A has moved the position of equilibrium a 
distance dr, the distance from the old position of equilibrium is a+dr. 


a i i i i i i a ee 
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The actual potential energy is the same, whether expressed as a power 
series with the old coefficients about the old position of equilibrium, or 
with new coefficients about the new position of equilibrium. Hence 


—Aa+a(x+dr)?+b(e+dr)3+e(xtdr)tt. .. 
=(a-+da)x?+ (b+db)x®+ (c-+de)et+. oak 


Since this is true for all values of x, we may equate coefficients of 
powers of x and obtain the equations 


da db | 
oe = 3b, WR CUCM, os ay ee) 
From eqns. (1) and (5), 
ldy 306 
ems 


This, with eqn. (4), gives 
_tdvlot) . 967, 
vdry oT 8a” 
If the distances between nearest neighbours in a crystal behave in the 
same way as the inter-atomic distance in this model, then 


be alin AAP O0F 
SERRE Sa8 5 é . . . . . . (6) 


where « is the volume coefficient of expansion and V is the volume of the 
crystal. 
Now d\inyv 
— Fe i ——Y . . . . . . . . (7) 


where y is the quantity used by Griineisen (1926) in his discussion of 
thermal expansion. It may be obtained experimentally from the 
relation 
is yw 8 
ca gen mae one ee” 
where K,, is the isothermal bulk modulus. If y is defined by eqn. (8), 
then eqn. (7) is not exact, but the approximation is reasonable in the 
present context. 
Equations (3), (6) and (7) now give 


5 4ca 


The dimensionless parameter ca/b? depends on the shape of the inter- 
atomic potential curve, but is not affected by stretching the curve 
uniformly parallel to either axis. It may be evaluated for any assumed 
type of inter-atomic force law. (Measurements of volume at high 
pressures give a and b, but are not in general accurate enough to give c.) 
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According to Slater (1939), a potential suitable for metals is that due 

to Morse : 
v=L exp {—2a,(r—rp)}—2L exp {—a,(r—") 5; 
where L and a, are constants, 7 is the distance between atoms and 1p 
is the equilibrium distance under zero external force. In this case ca|b? 
is 7/12, and 
8 
2Bk= ee 


A form of potential used by Mie (1903) and others is 
Mae 
ESP Capes =e mare 
ym yn 
where M, N, m and n are positive constants and nm is greater than m. 
In this case eqn. (9) gives 
2 mn— 2 
2B (5+ rata) 
Table 1 shows that the value of the quantity g=2Bk/y« is rather insensi- 
tive ton and m. Referring to eqn. (2), we therefore have finally, to the 
first order in 7’, 
Cy=k(1+qy«T') ee re se kd 


where q is expected to be not much less than unity. 


Table 1 


Morse 
potential 


0-71 0-70 ; : : 0-89 


Both y and « are slow functions of temperature, but here they have 
been regarded as independent of 7’, and the treatment has been to the 
first order in 7. Higher terms undoubtedly exist, and are likely to be 
of some importance near the melting point, but the limited purpose of 
the present treatment is to obtain an approximation for the linear term 
and compare it with the general trend of experimental results. 


§ 3. Discussion 


The thermodynamic relation between Cy and Cy, combined with 
eqn. (8), gives 
Cy=Cy(1+ yaT). 5. gy Sg eee em CH 


Equations (10) and (11) show that the increase of Cy above the Dulong 
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and Petit value is only a little less than the difference between Cp and 
Cy. The exact relation is of course dependent on the nature of the inter- 
atomic force law. To the first order in T, eqns. (10) and (11) combine 
to give 


1 /ac 
7 (Se) =O toe. eee ae 


In table 2, experimental values of the left-hand side for a few simple 
substances are compared with values of yx evaluated at twice the Debye 
temperature. 

The effect of using twice the Debye temperature rather than, say, room 
temperature is small but appreciable. The rather arbitrary choice 
could be avoided only by use of a higher order of approximation for Cy 
than that given by eqn. (10). For the metals listed, the electrons make a 
small contribution C. to the specific heat, and the second column of 
table 2 shows (0C./07')p/Cp calculated for one free electron per atom. 


Table 2. Units are 10-5 °c"! 


} 1 0Cp Vx Ratio 

(ar) z( aT ), ate Peo) Reference 
Pb 4 34 19 1-8 Douglas and Dever (1954) 
Cu 2 23 12 1-9 Landolt-Bornstein (1912-36) 
Na 4 85 35 2:4 | Dauphinee et al. (1955) 
NaCl — 43 24 1-8 Popov and Gal’chenko (1951) 
KCl — 39 21 1-9 | Skuratov and Lapushkin (1951) 
KBr — 38 22 1-7 Cooper (1953), Skuratovy and 


Lapushkin (1951) 


This has been subtracted from the observed (dC)/07')p/Cp to give the 
lattice contribution shown in the third column. The ratio in the fifth 
column is an experimental measure of 1-++-¢, which according to table 1 
is expected to lie between 1-7 and 1-9. . 

For sodium (0C'p/07)p/Cp is rather larger than suggested by this crude 
theory. Furthermore, it increases markedly near the melting point. 
The present discussion, having been taken only as far as terms linear 
in 7’, cannot indicate how much of this increase is to be expected in a 
perfect lattice and how much has to be assigned to other causes such as 
the thermal generation of crystal defects (Carpenter 1953). 

In recent discussions of the ‘linear chain’, Dugdale and MacDonald 
(1954) and MacDonald and Roy (1955) conclude that Cy should decrease 
with increasing temperature. They suggest that the difference between 
this behaviour and that of simple solids is a real difference between a 
linear chain and a solid body. This seems unlikely in view of the agree- 
ment obtained above with an even simpler model. | The relevant error 
in the general discussion of MacDonald and Roy is that they regard 
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their A, (equivalent to the a, b, c, etc. of the present discussion) as inde- 
pendent of volume, and therefore obtain too low a value for Cp. 


§ 4. CONCLUSION 


A crude theory predicts that above the Debye temperature Cy exceeds 
the Debye value by almost as much as Cy exceeds Cy, which roughly 
agrees with experimental results for several simple solids. 
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SUMMARY 
In a theory of thermal expansion above the Debye temperature, the 
increases of specific heat and Griineisen’s y with temperature are taken 
into account. The increase of expansion coefficient with temperature 
observed in several simple substances is shown to be about that expected 
for a perfect crystal. It is not evidence for the thermal generation of 
defects in the crystal. 


§1. INTRODUCTION 


THE volume coefficient of thermal expansion («) is normally observed 
to increase with increasing temperature, even well above the Debye 
temperature. This increase has been thought by some investigators 
(Lawson 1950, Uno 1951, Gertsriken 1954) to be a measure of the rate 
of increase of vacancy concentration with temperature. 

However, Griineisen’s theory of thermal expansion predicts an increase 
of « with increasing temperature for a perfect crystal, and has been 
used by Griineisen (1926), Nix and MacNair (1942) and Fischmeister 
(1956) to explain the increase observed for many simple substances. All 


‘these authors regard the two coustants of Griineisen’s equation (Q) and 


(m+n+3) in Griineisen’s notation) as independent constants to be 
determined empirically. They are in fact both related to the nature of 
the inter-atomic forces and hence to each other. When this relationship 
is taken into account it is found that Griineisen’s equation predicts only 


part of the observed increase of « at high temperatures. 


This discrepancy, noted by Fischmeister, was made clear by Fletcher 
(1957), who used the basic Griineisen hypotheses to obtain an explicit 
relation between«andtemperature. The increase of «at high temperature 
thus calculated for potassium chloride is about half the experimental 
rice theory, thermal expansion is related to the failure of 
Hooke’s law by noticing that the frequencies of atomic vibrations are 
dependent on volume. The failure of Hooke’s law is not, however, taken 
into account in dealing with the energy of the thermal vibrations, for 
which the Debye expression is used. This procedure is justified in 
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obtaining a value for « at low temperatures, but it will be shown that a 
better approximation is needed to obtain the variation of « above the 


Debye temperature. 
In the following treatment Griineisen’s law is obtained in a form which 


allows the theory of the variation of specific heat at high temperatures 
developed by Born and Brody (1921) and Schrédinger (1922) to be applied. 
An expression is thus derived for the variation of Griineisen’s y with 
temperature. The increase of thermal expansion with temperature is 
found to be about twice that given by Griineisen’s theory, and the result 
is compared with experiment for some simple crystals. 


§2. GRUNEISEN’s Law tn Exact Form 


Let p, V, 7’ and S be the pressure, volume, temperature and entropy of 
a fixed mass of a crystalline substance. Following Varley (1956) and 
Wilson (1941), we start with Maxwell’s relation 


ap\ _ (a 
OT }y \OV]m 
which may be written 


- (2) 4 eyes @ 
OV Ja\Ol/ oe NOT NOT 
Fg —— Cy 
Kpa2=— 7 (57), . . ° . ° . . (1) 


where « is the volume coefficient of expansion, K, is the isothermal 
bulk modulus and Cy, is the heat capacity at constant volume. With y 


defined by 
= —_ (F) ns ain | = . 
T\AV/ 5 On V/,’ 7) pes eels 


eqn. (1) may be written 


or 


yy 
—— VE, . . . . . . . . (3) 
or, using a well-known thermodynamic relation, 
bracers 
— VE. . . . . ° . . . (4) 


where K, is the adiabatic bulk modulus. Equations (3) and (4) are 
exact. To make them useful some information about y must be obtained 
from statistical mechanics. 


§3. THE VALUE OF y 
3.1. y in the Debye Approximation 
In the Debye approximation the entropy is of the form 
S=f(@/T) 
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where @ is a function of volume only. Hence, for constant S, 7’ is pro- 
portional to @, so that eqn. (2) becomes 


___aing : 


Thus y is independent of temperature in the Debye approximation. 

Since V and K’, are slow functions of temperature, eqn. (4) shows that: 
% i8 approximately proportional to C,. This is the original version of 
Griineisen’s law. 

When Griineisen’s formula for the temperature variation of K p is put 
into eqn. (3), together with y from eqn. (5), the resulting expression for « 
becomes identical, to the first order in 7', with the equation used by 
Griineisen, Nix and Fischmeister mentioned in §1. 


3.2. The Failure of Griineisen’s Equation 


The question is, whether the failure of this equation at high tempera- 
tures is caused by imperfections in the crystal or by a fault in the theory. 

It is not due to errors in the expression for K,,: according to eqn. (3), 
with y and C,, constant, « is inversely proportional to VK, but table 1 
shows that for several simple substances « actually increases appreciably 
faster than 1/VK >. 


Table 1 


1 (da ~ (SF) aC; (3) < 

ow mt), Ky\ oT }y Cy\ OT ] | Yo 

Pb 84 (Stokes and Wilson | 47 (Bridgman 1923) 15 4 is 
1941) ae He 

a rt 9.« 

Cu 47 (Nix and McNair 1941) Sone, on age bauney ul 2.2 


78 (Quimby and Siegel E 
Na 220 (Hagen 1883) tag : g 50 2-8 
NaCl 95 (Fischmeister 1956) 56 (Slater 1924) 19 2-0 
Kel 97 (Glover 1954) 38 (Slater 1924) 18 3:3 


Units are 10-5 °c-! Values of (dCy/0T),/Cy are computed at twice the 
Debye temperature. 


The cause of the failure must therefore be the prediction that yC, 
is constant at high temperatures. This prediction is not a result of the 
approximations peculiar to Debye’s theory, for, as may easily be shown, 
any theory which resolves the heat motion into that of a number of 
independent harmonic oscillators gives an expression for the entropy 
which leads to a constant value for y at temperatures high enough to 

i scillators fully. 
ar sie cigteanie ieietes the experimental results and the 
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predictions of Griineisen’s equation must be due either to imp s 
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in the crystal or to the use of the harmonic approximation. In the next 
section the latter possibility is investigated. 


3.3. First Approximation for the Temperature Dependence of yCy 


If the harmonic approximation were valid, (0C,,/dT)/Cy would, after 
subtraction of the electronic contribution, be small at high temperatures, 
but it is actually quite comparable with (0K 7/0T)/K 7. Since it is only in 
the harmonic approximation that y has been shown to be constant, we 
now try to assess the variation of y with temperature, using, in eqn. (2), an 
expression for the entropy which takes account of the anharmonic nature 
of the thermal motion. Born and Brody (1921) have shown that the 
Helmholtz free energy at high temperatures is, to a second approxima- 
tion, 

F=8RT (In (v/kT)—BkT}. 
Here v is the geometric mean of the frequencies of the normal modes at 
small amplitudes, and B is a quantity independent of 7’ which expresses 
the effect of departures from Hooke’s law. 

Since S=—(0F/0T),, it is easy to obtain 


Cy=T (eS/0T),=3R(1+2BkT) mee) Ai) 

and 
yOp=V(eOS/0V \p=38hyo(1--2BET e/y5) 9. ee ee 

where 
d In v 
and 
_adinB 

= inte te SUR emo eent Sa 


Thus yC; increases ¢/y,) times faster than C,, itself. Since B has not 
been calculated accurately, only a crude value of cis obtainable. Although 
« depends on the nature of the inter-atomic forces, which vary from 
substance to substance, it is possible to make some rough general 
estimates. 

In the theory of Born and Brody and Schrodinger, the energy is stated 
in terms of the normal momentum and position coordinates, p, and q, of 
the crystal in the form 


i> (p+ w?q77) +f5(9;) +ha(9:)- 


The w, are the angular frequencies of oscillation of the normal modes at 
infinitesimally small amplitudes. The quantity fs is the sum over all 
values of 2, j,k, of terms b,,, 9; q; 7, where the b,,,, are unknown coefficients 
which depend on the nature of the inter-atomic forces. Similarly f, 
is the sum over #, 7, kh, 1, of ¢;5419;4;0;¢;. It may be shown that the value 
of B is then 


Bay —2P inn Casket 


A , > 
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where the first summation includes only those terms for which the 
denominator can be expressed in the form w ,?w,*w,2 and the second 
summation includes only those terms for which the denominator has the 
form w ,2w,?. 

If we assume that all the coefficients of the second-, third-, and fourth- 
order terms in the expression for the potential energy depend on volume 
in the same way, then 

dinB dino? > din 
ee “ei v= dae 
Hence, from eqns. (8) and (9), 


e/yo=2. ie hire shee fH SCLO) 
The derivation of this equation takes account of the coupling between 
atomic vibrations, but makes the crude assumption that all the w,?, 
b,;, and ¢;;,,; alter in proportion when the volume is altered. Another 
approach is to ignore the coupling and treat a two-atom model, using an 
empirical form of inter-atomic potential to determine the volume 
dependence of the coefficients. It has been shown in a previous paper 
(Eastabrook 1957) that the two-atom model allows a reasonable estimate 
of the quantity B. From eqns. (1), (3) and (5) of that paper one easily 
obtains, the same notation, 


ely 15 Boe ons 15 3ca 
Wyn BAI @.3b? 1259) \ 16 4b?) 


Table 2 shows values of ¢/y,) obtained with a Morse potential, and with a 
Mie potential using various values of the exponents m and n. It will be 


Table 2 
n iS; 9) 5 9 9 Morse 
m 1 1 3 3 6 potential 
</Vo 2-48 2-44 2-86 2°86 3:16 Seal 7 


noted that ¢/y» is not very sensitive to m and x. Thus, from both table 
2 and eqn. (10) we expect </y) to be about two or three ; that is, yC, 
should increase two or three times faster than C,,.. A comparison (table 1) 
of the rates of increase of (VK)! and C,, shows that according to eqn. (3), 
about half the increase of « should be related to the increase of yC, 
and the rest to the decrease of VA‘ p. 

From eqn. (3), 


1 (eer) =+(5r) ! — (Sar) rer 
TOg\ Ol je ON) Kn \ OLY jr 


Thus an experimental measure of </yo is obtained by comparing the 
quantity on the right with (00,,/0T),,/Cy after ailowing, where necessary, 
for the small electronic contribution. Independent measurements of the 
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quantities needed in eqn. (11) are rarely concordant, and the figures for 
</y, in table 1 should be regarded merely as evidence of a general trend. 
They are, as predicted above, about two or three. 


§4. CONCLUSION 


The variation of y and C,, with temperature must be taken into account 
when applying Griineisen’s law to the increase of thermal expansion at 
high temperatures. A rough calculation of the variation of yCy shows 
that the increase of « observed in each of several simple substances is 
such as would be expected for a perfect crystal. Thus the increase of 
thermal expansion at high temperatures cannot be regarded as evidence 
for the thermal generation of defects in the crystal. 
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SUMMARY 


The internal friction of cold-worked polycrystalline copper has been 
measured in the temperature range 20° to 300°K, and a study made of 
the maximum in the curve of internal friction versus temperature first 
observed by Bordoni. The measurements indicate how the internal 
friction varies as a function of the amount of cold work, and also the 
effect of impurities in the metal and of neutron irradiation. The results 
are discussed with particular reference to recent theoretical treatments 
by Seeger. 


§ 1. INTRODUCTION 


THE internal friction of several cold-worked metals in the temperature 
range from 4° to 300°K has been measured by Bordoni (1949, 1954), 
who showed the existence of a peak in the curve of damping against 
temperature. Bordoni’s principal results, obtained at frequencies 
between 10 and 40 ke/s, are shown in fig. 1, the internal friction being 
expressed in terms of the quantity 1/Q (the logarithmic decrement 6 
divided by z). The largest peak occurs for copper at a temperature of 
about 90°x. For aluminium, lead and silver the damping is smaller ; 
the peaks are less pronounced, and occur at somewhat different tempera- 
tures. The existence of such peaks has subsequently been confirmed by 
the present authors (1955) in the case of copper, and by Bommel (1955) 
and Hutchison and Filmer (1956) who made measurements at frequencies 
in the megacycle range on single crystals of lead and on polycrystalline 
aluminium respectively. ; 

As no completely satisfactory explanation of these peaks was available, 
the present work was undertaken in order to obtain more systematic 
information ; in particular to investigate the internal friction. as a function 
of cold work more closely than had been done previously, All the 
present experiments have been made on copper, as this metal appeared 
to show the largest effect Measurements have also been made to 
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investigate the effects of neutron irradiation and of the introduction of 
impurity atoms. Two short notes giving a brief account of some of our 
results have already been published (Niblett and Wilks 1955, 1956). 


Fig. 1 
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Internal friction in metals (Bordoni 1954), 


A. Chemically pure polycrystalline copper. (a) machined ; (b) after 10 hours 
annealing at 150°c. 

B. 99-9% pure aluminium. (a) after 1 hour anneal at 140°c; (b) after 2% 
permanent strain. 

C. Lead. (a) chemically pure; (6) commercial, 

D. Silver. (a) machined ; (b) after 14 hours anneal at 200°c. 


§2. APPARATUS 


The principal features of our experimental arrangement are illustrated 
in fig. 2. The specimen, a rectangular bar 8 cm long, 1 cm wide and 
about 1-8mm thick, is supported on the edges of two razor blades, 
which are positioned at the nodes of fundamental flexural vibrations. 
The bar is driven at its resonant frequency by applying an alternating 
voltage to an electrode placed beneath its centre, the separation of the 
electrode from the specimen being about 0-3mm. To detect the resulting 
motion of the bar, two electrodes are placed one at each end of the 
specimen, thus forming small parallel-plate condensers, each with a 
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capacity of about 3uur. The movement of the specimen produces 
changes in the capacity of these condensers, whose magnitude is propor- 
tional to the amplitude of the vibrations. These variations in capacity 
are measured by a circuit in which the condensers are connected in 
parallel with the tank circuit of a quartz oscillator ; motion of the speci- 
men produces a modulation of the output of the oscillator, and this is 
displayed on a cathode-ray oscilloscope. A measurement is made by 
observing the free decay of the oscillations on switching off the drive 
voltage, either by photographing the oscilloscope trace or with a level 
recorder. (A level recorder with a suitable logarithmic potentiometer 
converts an exponentially-decaying signal into a linear pattern on a 
moving strip of paper, and the internal friction is given immediately 
from the slope of the trace.) 


Fig. 2 


Dewar 
vessel 


Platinum 
Resistance 
hermometer 


Apparatus for measuring internal friction at low temperatures. 


The absolute value of the strain amplitude may in principle be 
calculated from the driving voltage and the separation of the plates, but 
this is rather difficult as the gap is small and not always quite uniform, 
Also as the strain varies throughout the specimen, the experiments 
measure the internal friction for a range of strain amplitude ; we therefore 
indicate the amplitude of the vibrations by giving an estimate of the 
mean strain. In spite of considerable uncertainty in the absolute value 
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of the strain, it is of course possible to obtain accurate values of the 
relative amplitudes on a given specimen in situ. 

The razor blades which support the specimen are clamped rigidly 
in small brass blocks, which are screwed into a groove in a brass base- 
plate ; this is attached by four pillars to a further brass block B, weighing 
about a kilogram, suspended by three wires from the top of the cryostat. 
(The fine copper leads soldered to the wires at X ensure that any heat 
coming down the wires is taken up by the liquid in the dewar.) The 
cryostat itself is in turn mounted on a massive iron framework embedded 
in a block of concrete, the whole assembly being designed to reduce the 
effect of external vibrations to a minimum. The specimen and brass 
block, B, are enclosed by a vacuum jacket, J, which in turn is surrounded 
by a second vacuum jacket immersed in liquid hydrogen or oxygen. 
Heating coils are attached to the brass block, B, and to the outside of the 
jacket, J, in order that measurements can be made over the whole tempera- 
ture range from 20°K to 300°K. Helium exchange gas at a pressure of 
about 2mm of mercury is introduced into the inner vacuum space to 
establish thermal equilibrium within the jacket J; this space is, however, 
evacuated before a measurement is made, in order to eliminate damping 
due to the gas. Temperatures are measured by means of a platinum 
resistance thermometer, wound on a brass former attached to the under- 
side of the brass block. 


§ 3. RESULTS 
3.1. Preliminary 


We have subtracted from all our results the damping due to the thermo- 
elastic effect, using an expression due to Zener (1937) which has been 
verified experimentally by Bennewitz and Roétger (1938) for german silver 
reeds, by Zener et al. (1938) for copper, and by Berry (1955) for «-brass. 
In all cases the correction is a small fraction of the total damping, being 
less than 2% in the region of the peaks. Unless otherwise stated the 
measurements were made on high purity oxygen-free polycrystalline 
copper (99-999°%,) supplied by the British Non-Ferrous Metals Research 
Association. All the specimens were first annealed in argon for one 
hour at 600°C, and then plastically deformed at room temperature by 
longitudinal extension in a Hounsfield tensometer. The ends of the 
specimen, which had been gripped in the jaws of the tensometer, were 
cut off before the internal friction was measured. In a preliminary 
experiment a specimen was strained 4%, and its internal friction 
subsequently measured at room temperature as a function of time. 
A large fraction of the increase in damping produced by the deformation 
disappeared within a few hours; this is the so-called Koster effect first 
observed in brass and steel (Forster and Koster 1937, Késter and Rosen- 
thal 1938). All the specimens in the present experiments were therefore 
left at room temperature for several days after deformation before 
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measurements were made, so that anv f: i 
: , y further changes in the damping, 
associated with the Koster effect, would be negligible. oa 


3.2. The Effect of Cold Work 


The temperature dependence of the damping of one annealed specimen 
prior to deformation is shown in fig. 3. The two curves were obtained 
during separate experiments on the same specimen at strain amplitudes 
in the ratio of 2 to 5, the value of the strain amplitude being of the order 
of 10- *. At higher values of the strain amplitude the damping increases 
rapidly with increasing amplitude at all temperatures ; hysteresis effects 
are observed and it becomes difficult to obtain reproducible readings 
(Niblett 1956). | 


Fig. 3 


0 50 100 150 200 
Temperature °K 


The internal friction of annealed copper. Strain amplitudes in 
arbitrary units : + 2, O5. 


The effect of varying amounts of cold work up to 8-4% extension is 
shown in fig. 4; all the measurements were made at a strain amplitude 
of the order of 10-7, and at a frequency of approximately 1100 cycles. 
per second. The most obvious feature of the curves is the existence 
of two peaks in the internal friction at temperatures of about 32°K and 
75°K which become more pronounced in the more strongly cold-worked 
specimens. For amounts of cold work greater than 2% the peak at the 
higher temperature is similar in shape to that observed by Bordoni (1949, 
1954) for copper, although his measurements do not seem to have been 
sufficiently detailed to resolve the second peak. (However, Bordoni’s 
results do show some irregularity in the region where this peak might 
be expected.) The shape of the curve for copper strained 0-5°% is similar 
to that obtained by Bordoni for silver, while that for copper strained 
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Fig. 4 
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The internal friction of copper. 


A. Strained 0-1%. B. Strained 0-5%. 
C. Strained 2-:2%. D. Strained 8-4%. 


Fig. 5 


Temperature °K 


Amplitude dependence of the internal friction of copper, strained 2-2%%. 
Strain amplitudes in arbitrary units: 41, + 10, ©100. 
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0-1% bears some resemblance to Bordoni’s curve for lead. Hence the 
marked variation in the shape of Bordoni’s curves may arise not so 
much from the difference in materials but rather from the metals being 
subjected to different amounts of cold work. 

In the more strongly cold-worked specimens the damping is to a 
first approximation independent of the strain amplitude ; this is illus- 
trated in fig. 5, which shows the internal friction of a specimen strained 
2-2°%, measured at three strain amplitudes in the ratio 1: 10: 100 (the 
absolute values being in the range 10-®-10~6). This absence of strong 
amplitude dependence, together with the fact that the stresses applied 
during a measurement have virtually no effect on the state of strongly 
cold-worked specimens, results in these measurements being much more 
reproducible than those on annealed copper. 

At first, both peaks increase in amplitude (at a similar rate) with 
increasing deformation, and appear to be superposed on a background 
of friction arising from some different mechanism, presumably that 
responsible for the damping in the annealed state. Ultimately at exten- 
sions of more than about 2°, the amplitude of the peaks becomes approxi- 
mately independent of the amount of pre-strain. It is to be noted that 
the peaks occur at approximately the same temperatures for each speci- 
men. It seems that the friction of specimens which have suffered only 
small pre-strains may be resolved into two components consisting of 
peak and background, and in this case also the shape and position of the 
peaks appear similar to that at the higher deformations. Therefore we 
propose to treat the mechanism responsible for this peak as a separate 
process, apart from that responsible for the background, and it is with 
this mechanism that the present paper is principally concerned. The 
question of the background friction is considerably more complex and 
is still not understood in detail even at room temperature. 


3.3. Effect of Frequency - 

All the measurements so far described were made at a frequency of 
about 1100 cycles per second. Although it is of interest to obtain 
measurements at different frequencies, it is difficult to vary the fre- 
quency over a wide range with the present apparatus because this involves 
changing the thickness of the specimen. At higher frequencies the 
specimen rapidly becomes too thick and stiff to drive; at lower fre- 
quencies it becomes too thin to handle. Measurements have, however, 
been made on a thinner bar with a resonant frequency of 380 cycles per 
second, which had been strained 9-1"%. The friction was similar to that 
for similarly strained thicker specimens, except that the internal friction 
peaks now occurred at about 67°K and 30°K (Bordoni’s measurements at 
30 ke/s show a peak at about 90°K.) . : 

Thus in the more heavily cold-worked specimens, the damping is 
approximately independent of strain amplitude, the maximum es 
is substantially independent of frequency, and the temperature 7’ where 
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the peak occurs varies with frequency. This suggests that the internal 
friction may arise from a relaxation mechanism, and one might expect to 
find a relation of the form 

f=A exp (—W/kT) Jt ee Lasts 
where f is the frequency of measurement, A a constant, and W the 
activation energy for the process (see for example, Zener 1948). Hence, 
a plot of the logarithm of the frequency against the inverse of the tempera- 
ture at which the peak occurs should give a straight line whose slope is 
the activation energy. The available information consists of our results 
at about 1100 c/s and 380 c/s, and Bordoni’s measurements at about 
30 ke/s. These do not lie very closely on one line, possibly because the 
peak cannot be attributed to one single activation energy, as is also 
suggested by the facts that the peaks are broader than might be expected. 
However, taking the best line possible one obtains a value for the activa- 
tion energy of about 0-1 ev. 

A further criterion for relaxation processes is that there is a change 
in the elastic modulus, #, on passing through the region of the peak, 
such that 4H#/H=1/Q,, where Q,, is the value of @ at the peak. There 
is some uncertainty in obtaining exact values of the modulus change 
associated with the peaks because of the presence of the background 
friction, but a corresponding modulus change of this order is observed 
in the more heavily cold-worked specimens. 


3.4. Effect of Annealing 

The effect of annealing at various temperatures was observed using 
specimens which had been strained 0:1°%% and 8-4%, and therefore 
exhibited damping typical of the background and the peaks respectively. 
All the annealing was carried out in an atmosphere of argon. An anneal 
for one hour at 200°c had practically no effect on the internal friction 
of the specimen strained 0-1%, but a subsequent anneal for one hour at 
250°C reduced the damping over the whole temperature range to a value 
as low as that for an unstrained specimen. In the case of the specimen 
strained 8-4°, an anneal for one hour at 180°c reduced the heights of the 
peaks by about 40°%, without changing the general shape of the curve 
(fig. 6) ; this result is similar to that obtained by Bordoni after an anneal 
at 150°c for 10 hours (fig. 1 (@)). A subsequent anneal for one hour at 
250°C had little additional effect on the damping. After a further 
anneal for one hour at 350°c, the peaks had disappeared and the curve 
of damping against temperature was similar to that for a fully annealed 
or slightly strained specimen. It is also to be noted that the measure- 
ments shown in fig. 5 were obtained with the same specimen as that used 
for fig. 4 (c) after a lapse of several months ; thus no appreciable annealing 
occurs at room temperature. Hence it appears that the peaks are not 
to be associated with the motion of free defects created by the deformation. 
This suggests that the relaxation mechanism responsible for the internal 
friction may be the motion of segments of dislocation line, with a com- 
paratively low activation energy, of the order of 0-1 ev. 
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3.5. Effect of Impurities and Neutron Irradiation 


Measurements have been made on copper containing 0-0026°, bismuth. 
and 0-032°%, phosphorus, and the results for specimens strained OLE: 
and 5-5°, are shown in fig. 7. The curves are rather similar in shape 
to those for correspondingly strained pure copper, but in each case the 
damping at all temperatures is reduced by a factor of about ten. A. 


Fig. 6 
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Effect of annealing on the internal friction of copper. O After straining 8-4. 
+ After 1 hour anneal at 180°c. [After 1 hour anneal at 350°c. 


Fig. 7 
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The internal friction of copper, containing 0-0026% bismuth and 0-03 
+ Strained 0:17%. © Strained 5-5%,. 
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similar reduction in the friction has been observed after neutron irradia- 
tion as is shown in fig. 8 for a specimen which was first strained 6-6°% and 
then irradiated to a flux of 5X 1018 n.v.t. in the pile at the Atomic Energy 
Research Establishment, Harwell. There appears to be a slight shift 
of the temperature of the peak in these experiments, of about 5° to lower 
temperatures, but to a first approximation the position and shape remain 
unchanged. Thus even though the motion of the dislocation lines is 
considerably reduced the activation energy for the process remains 
substantially unchanged; and this suggests that it does not depend 
markedly on the distance between pinning points. Also the effect of 
cold work is to produce a more complex dislocation net with shorter 
lengths of free lines, yet the activation energy remains unchanged. Thus 
the length of dislocation lines responsible for the internal friction peak 
cannot be specified either by the distance between the nodes of a disloca- 
tion net or between pinning points associated with impurity atoms or 
points defects. (This conclusion is at variance with an explanation 
of these peaks offered by Mason 1955.) Clearly some form of unit 
displacement of the dislocation line is involved ; and there appear to be 
two possibilities. Either a bulge is formed by thermal activation as 


Fig. 8 
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Temperature °K 
The internal friction of copper, strained 6-69 and subsequently neutron 
irradiated (5 x 10!8 n.v.t.) 


discussed in $4, or the line becomes pinned and unpinned from some 
point defects. In the latter case the magnitude of the friction would 
depend critically on the amount of impurity present. Nearly all our 
work was done with 99-999, high purity copper in which one might have 
supposed that the impurity content was quite low. Another experiment 
was performed on copper from a different source with a significantly 
better specification ; the specimen was annealed in vacuo, strained 5-3, 


and had a resona V ' ae 
esonant frequency of about 750 c/s. Apart from the presence 


of a third peak at about 200°k which we refer to later, the friction is 
similar to that observed before, although the height of the peak is greater 
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(fig. 9). Bordoni’s results at 30 ke/s also give the same magnitude of 
peak for nominally pure copper from some quite different source. If 
unpinning from impurities were responsible for the friction, it is rather 
remarkable that all these experiments give such similar results. On the 
other hand, such a situation is to be expected on a theory due to Seeger 
which is described below. 

Fig. 9 
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The internal friction of the high purity copper, strained 5-3°%. 


§ 4, DiscussIon 
4.1. Seeger’s Treatment 


According to Seeger (1955 a), internal friction will arise from the motion 
of dislocation lines lying parallel to lattice directions. Consider a length 
of dislocation (AB in fig. 10 (a)) lying in a position of minimum energy 


Fig. 10 


Seeger’s dislocation mechanism. 
(a) Dislocation in position of minimum energy. 
(b) Bulge in dislocation. (Full straight lines represent positions of 
minimum energy, broken lines.represent positions of maximum energy.) 
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along a close-packed direction in a lattice, and suppose that an applied 
stress is tending to move it forward in a direction parallel to itself. lim 
the dislocation acts as a rigid rod, a stress equal to the Peierls’ force 
will be required to move it over the potential barrier to the next equili- 
brium position CD, but if the line is flexible it may be able to move 
forward much more easily. Thus, if a small bulge (XY in fig. 10 (b)) is. 
formed as the result of thermal fluctuation, a quite small stress will cause 
it to spread out sideways, as the potential barriers opposing motion in 
directions parallel to AB and CD are very small. The length of a bulge 
formed by thermal activation will be determined by the fact that its. 
energy of formation is a minimum; a longer bulge would require more 
energy to surmount the Peierls’ barrier, while a shorter one would involve 
a greater increase in line length and energy. A bulge of this type may 
be regarded as a pair of kinks of the kind which arise when a dislocation 
passes from one lattice row to the next. In the absence of any stress. 
the two kinks will attract and annihilate each other, but an applied 
stress will tend to move them apart and increase the slipped area. For 
a given applied force there will be a certain critical separation of the 
kinks derit above which they will separate further, and below which they 
will come together. Internal friction arises under the influence of an 
applied stress because of the formation of these bulges; the motion of 
the two constituent kinks is approximately reversible for small displace- 
ments, but as soon as the critical distance derit is exceeded they will move 
apart rapidly and irreversibly. The internal friction will be of the general 
form associated with relaxation processes, as for example is the internal 
friction arising from the motion of solute atoms in solid solutions sub- 
jected to stress induced ordering. 

The activation energy for the relaxation process responsible for the 
internal friction will be W the energy of formation of a bulge. In his. 
initial treatment Seeger (1955 a) took W to be 2w where w is the additional 
energy associated with a single kink in a dislocation which is otherwise 
parallel to a lattice direction. The value of w is determined by specifying 
the form of the dislocation by the differential equation 


dy 2ry 
isa =bo, (1 cos =a 


where 18 is the energy of unit length of line, a the lattice spacing, b the 
Burgers’ vector, and a, the Peierls’ stress ; hence Seeger (1956) obtains. 


the value 
O== 2a JCS) + 
7 aT 


However, this treatment is admittedly only an approximation as it 
ignores the precise form of the energy barrier which must belsurmounted 


} This expression which appears to be th : 
from that previously given by Read (1953) e correct one differs by a factor 2° 
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‘during the formation of a bulge. A rather more extensive but. still 
-approximate treatment (Seeger 1956) leads to the expression 


W=v{1-+} log (160,/70)] 


where o is the applied stress: this implies that a factor of 100 on the 
magnitude of the amplitude changes W by about 20° and so shifts the 
position of the peak by about 15°. As no such shift is observed (fig. 5), 
it appears that the approximation is not a very good one. However 
Seeger et al. (1957) have recently reported the results of a much more 
complete treatment in which the relaxation time is derived from the 
theory of stochastic processes. This treatment, which has not yet been 
published in detail, is rather involved so that the expressions for the 
activation energy and frequency term in eqn (1) cannot be given in an 
explicit form ; reference must therefore be made to the aforesaid paper. 
In the next paragraph we compare this treatment with our experimental 
results. 

The most satisfying aspect of Seeger’s treatment is that his mechanism 
accounts for a relaxation process with a quite closely defined activation 
energy, and in the latest treatment this depends only very slightly on the 
strain amplitude for strains of the order used in the experiments. The 
treatment also accounts for the presence of two distinct peaks in the 
internal friction—temperature curves. Only dislocations lying parallel to 
the principal lattice directions will move in the manner discussed ;  dis- 
locations lying in other directions may move easily under the action of 
much smaller stresses (Mott and Nabarro 1948) and the relaxation 
mechanism will not arise. As dislocations lying parallel to lattice 
directions in a close-packed f.c.c. lattice may have Burgers’ vectors 
-corresponding to either a screw or mixed screw-edge orientation, we may 
expect two activation energies, one for each orientation. 

The initial increase of the internal friction with cold work is to be 
associated with the increase in the number of dislocations, and there are 
at least two possibilities as to why the friction does not continue to 
increase for pre-strains greater than 2%. Possibly the density of dis- 
locations continues to increase, but the effect of their interactions becomes 
‘so pronounced that they are less liable to lie along lattice directions 
(Seeger 1956). Alternatively, the number of dislocations may increase 
much more slowly at the higher stresses, the increased strain arising 
primarily to the motion of dislocations already present surmounting 
barriers at which they were previously held up. In either case the 
behaviour of the internal friction should be reflected in the stress-strain 
curve, and Feltham and Meakin (1957) have observed a change from 
linear to parabolic hardening in polycrystalline copper for a strain of 
about 14%. (This result being substantially independent of grain size.) 

In order to calculate the value of the activation energy from this theory, 
it is necessary to assume expressions for the energy (Hy) and the inertial 
mass (m) of unit length of dislocation, and for the Peierls’ stress o,. As 
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the value of this last parameter is still rather uncertain compared. with 
the other two, Seeger et al. assume that Hy=}Gb* and m=E,/v2 (v; is the 
velocity of shear waves in the metal), and use the experimental value of 
the activation energy to calculate a value for o,. For the strain ampli- 
tudes used in the present experiments, the results of such a calculation 
are substantially independent of the amplitude, and give o,=2-5 x 1034G: 
If we wish to throw the experimental results into the form of eqn. (1), 
then Seeger’s treatment shows that the term A is weakly temperature 
dependent but that to a reasonable approximation is about 101°. This is 
considerably lower than the value obtained in previous treatments which 
regarded the dislocation as a rigid rod oscillating in the potential trough, 
but is still appreciably higher than the experimentally determined value 
of about 10°. The magnitude of the maximum friction depends on the 
total length of dislocation in the crystal, on the fraction which is lying 
parallel to lattice directions and on the length of individual runs of dis- 
locations. Although it is hardly possible to make an exact estimate, it 
seems that if the total length of line in cold-worked material is taken to 
be 101° cm~2, then only about 0-1°% of the dislocations would contribute 
to the damping. This would appear to be a reasonable figure in view of 
the small size of the Peierls’ force compared with the other forces acting 
on the dislocations. It may be noted that electron micrographs (Hirsch 
et al. 1956) do not show much evidence of dislocations lying parallel to 
lattice directions but it must be remembered that the scale of the 
mechanism proposed by Seeger is quite small. Thus the critical length of a 
bulge is estimated to be about 60 atoms, so that even if the mechanism 
only occurs when runs considerably longer than this are parallel to lattice 
directions, lengths of a few hundred 4 would be sufficient to produce the 
effect. These could well be present and not show up clearly in the 
existing photographs. 


4.2. The Effect of Impurities and Neutron Irradiation 


The internal friction is much reduced at all temperatures in copper 
containing quite small quantities of bismuth and phosphorus. Such an 
effect is well established at room temperatures (e.g. Marx and Koehler 
1950), and must arise from pinning of the dislocations by the impurities. 
Although a similar explanation will still be valid at lower temperatures in 
the region of the peaks, two points call for comment. If, as seems 
probable, the friction at room temperature arises from the motion of 
dislocation segments about 1 long, then an impurity concentration 
of 1 in 10* may have a considerable effect on their mobility. At first 
sight, such a concentration appears hardly adequate to reduce friction 
arising from the motion of segments only 60 atoms long. However 
during the cold work the dislocations move through the specimen and 
will interact with any impurity atoms they encounter; thus the con- 
centration of impurities along a dislocation will be greater than the 
average throughout the specimen. Secondly, if we accept the above 
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explanation of the two peaks, the results imply that pure screws as well 
as mixed dislocations may be pinned by impurity atoms. This may 
arise either through their splitting into partials, or due to electronic 
interaction (Stehle and Seeger 1956). The higher damping observed in 
the nominally purest copper specimen (fig. 9) may of course be explained 
by a reduction in the number of pinning points. 

The effect of neutron irradiation on a specimen is to reduce the friction 
of deformed specimens to values rather similar to those observed in impure 
copper, and it appears that the dislocations are now pinned by inter- 
actions with point defects created by the irradiation. Measurements of the 
change in electrical resistance produced by irradiation (e.g. Marx etal. 1952) 
together with an estimate of the resistivity associated with one vacancy 
or interstitial suggests that the number of point defects introduced will 
be roughly of the same order as the concentration of impurity atoms in 
the specimens previously discussed. The pinning is presumably due to 
jogs which result from interactions between dislocations and point defects. 


Fig. 11 
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The internal friction of copper, neutron irradiated (5x 10'% n.v.t.) and sub- 
sequently strained 5-1%. 


and which restrict the motion of the dislocations (Seeger 1955 b).  Irradi- 
ation prior to deformation does not reduce friction, but this is hardly to 
be expected because most of the dislocations responsible for the friction 
are only introduced by the cold work. In fact the friction has a value 
somewhat higher than that of the unirradiated copper (fig. 11), and it may 
be that this is a secondary effect of the increased hardness brought about 
by irradiation (Billington and Siegel 1950). This is accompanied by a 
much coarser mode of slip (Jamison and Blewitt 1952), and it may be that 
for a given amount of slip there is less interaction between the dislocations 
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themselves. It is also to be noted that the specimen deformed after 
irradiation shows a small broad peak in the region of 200°K ; this peak 
which appears common to all specimens so treated is similar to that 
observed in the nominally purer copper. 


4.3. The Peierls’ Stress 


It would appear from the above discussion that Seeger’s mechanism 
gives a good qualitative account of the experimental results, but quanti- 
tatively does not lead to the correct magnitude for either the activation 
energy or the frequency factor. It is true that the value for the Peierls’ 
force is similar to ones calculated theoretically (Dietze 1952, Seeger 1955 a) 
but as is well known this latter figure depends very critically on the width 
of the dislocation (e.g. Foreman et al. 1951). The latest and most complete 
treatment (Seeger ef al. 1957) leads to a value for the Peierls’ stress con- 
siderably lower than those obtained previously but which is still consider- 
ably higher than the yield stress of single crystals (for example, Blewitt 
-et al. (1954) quote values of about 10~4@). 

Although the Peierls’ stress may be related in some way to the yield 
stress of single crystals it does not equal it, as the Peierls’ stress is the force 
needed to move a length of dislocation line only if the line lies parallel to 
a close packed direction and if the dislocation moves as a rigid rod. 
Our particular concern is whether the Peierls’ stress can be higher than 
the yield stress in material where the yield stress depends neither on 
impurities nor interactions between dislocations. To obtain appreciable 
‘plastic flow the stress must be at least sufficient to operate Frank—Read 
sources, and the critical stress is usually calculated from the line tension 
of the dislocation, the effect of potential wells associated with the lattice 
lines being ignored. This last effect is small for dislocations which run 
across lattice directions (Mott and Nabarro 1948) but in general some 
segments of a dislocation will lie parallel to lattice directions. If these 
segments are to move forward, and not hold back the rest of the line, 
then it might appear that a stress equal to the Peierls’ stress must be 
applied. (Dislocations which are already present in a crystal and which 
he across lattice directions will of course be able to sweep out a 
considerable area before being held up at a Peierls’ barrier, but the total 
strain which can be produced in this way is negligible.) However there 
are at least two possibilities whereby the dislocations might cross the 
barrier under a lower applied stress than the Peierls’ stress. The 
mechanism discussed in § 4.1 would enable the barriers to be crossed at 
high temperatures by the formation of bulges, but would imply a much 
stronger dependence of the yield stress on temperature than is observed. 
Nevertheless according to the figures given by Seeger et al., the frequency 
at which bulges are formed increases considerably for values of a/o 
greater than about 1°%, so that even at low temperatures the dislocation 
might be able to surmount the barriers comparatively easily. Another 
possibility (Granato, private communication) arises from the fact that 
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a dislocation not originally parallel to lattice directions will have moved 
some distance before having to surmount a barrier, and may have 
sufficient kinetic energy to carry itself over. 


§ 5. CONCLUSION 


It has been shown above that the mechanism proposed by Seeger 
accounts qualitatively for our results, but leads to rather high values 
for the activation energy and the jump frequency. It may be that a more 
refined calculation, say taking into account the fact that the dislocations 
will be split into partials, would lead to better agreement, but one should 
also consider if there is any alternative mechanism which could be 
responsible. Our results seem to indicate that the friction arises from 
the motion of lengths of dislocation line which can move with a small 
activation energy, but only over small distances. As mentioned pre- 
viously it is possible that the process responsible consists of dislocations 
being pinned and unpinned from impurity atoms, and such a mechanism 
may well have an activation energy of about 0-lev. Although this is in 
some ways an attractive suggestion, the experimental evidence is rather 
against it. Thus the magnitude of the friction would presumably be 
very dependent on the impurity concentration, but our specimens of 
copper from two sources and those of Bordoni all give the same order 
of friction. Also there is no marked amplitude dependence of the friction 
as might be expected if it arose from unpinning. Thus the evidence 
suggests that a mechanism of the type proposed by Seeger is responsible 
for the two peaks, and this could probably be confirmed by measurements 
on zone annealed copper when it becomes available. Seeger’s theory 
does not explain the peak which is observed in some specimens at about 
200°K, but the fact that it is not always observed (in contrast to the 
Bordoni peaks) implies that it arises from some other mechanism. 

Finally it is worth noting that a very marked feature of the Bordoni 
peaks is their reproducibility and independence of strain amplitude, both 
of which characteristics are in marked contrast to the behaviour of the 
background friction. In this respect polycrystalline material is markedly 
superior to single crystals ; although both peaks have been observed in 
single crystals (Niblett 1956) they are superposed on a greater amount of 
background friction and their analysis is correspondingly more difficult. 
Thus the present relaxation process, once fully understood, would appear 
to offer a useful way of observing the motion of dislocation lines under 
well defined conditions. For example, it is possible that the effect of 
neutron irradiation on internal friction might be more profitably studied 
by its effect on these peaks rather than by using annealed material. 
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SUMMARY 


X-ray and density studies have been made on single crystals of 
quartz irradiated with fast pile neutrons at approximately 100°. 
The structural effects of the irradiation follow at least two processes 
with doses up to 1-2 102°neutrons/em?. Crystals irradiated with 
less than 3101! neutrons/em? expand anisotropically in a manner 
that is similar to the thermal expansion of unirradiated quartz. In 
this dosage range the volume increase is the same, as determined by 
x-ray and hydrostatic methods, and little or no distortion is observed 
in x-ray diffraction patterns. The model which seems to fit these results 
best consists of point defects and slightly disordered regions which result 
in an elastic distention of the lattice. 

After doses in excess of 3 101° neutrons/cm? the damaging process is 
more complex since the volume change indicated by x-ray measurements 
of the host lattice becomes larger than the bulk volume change determined 
hydrostatically. After irradiation with 8 x 1019 neutrons/em? an in- 
homogeneous shear strain is observed in the extreme skewing of the 
(22-0) reflection which is interpreted as being due to the crowding of 
oxygen ion interstitials into the open c-axis channels. Amorphous 
regions enmeshed in the host lattice are identified in crystals that still 
possess long range order after irradiation with 7x 10! neutrons/cm?. 
The diffuse x-ray scattering halo associated with the amorphous material 
changes peak position, shape, and intensity with increasing dosages of 
irradiation. It is suggested that these changes result from a reorientation 
of disordered regions which is assisted by the increased expansion of the 
lattice and the presence of ruptured silicon—oxygen bonds. A second 
kind of diffuse scattering appears at these high doses which cannot be 
accounted for on the basis of a thermal origin but is probably related to 
scattering resulting from inhomogeneous lattice strains. 

In a complementary study it was found that coesite, a high density 
crystalline silica remained stable under neutron irradiation dosages which 
completely disorder quartz. These results indicate that the rate of 
damage in silica solids is structure dependent. 


+ Communicated by the Author. 
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§ 1. INTRODUCTION 


Tu influence of different types of lattice imperfections upon the low 
temperature thermal conductivity of solids (Klemens 1951) has provided 
a sensitive method for determining the nature of disorder in neutron 
irradiated quartz (Berman 1951). As a result of these earlier investiga- 
tions, the presence of defects with atomic dimensions and more complex 
regions of disorder were first indicated, and Berman’s study also revealed 
that the low temperature thermal conductivity curve for neutron irradiated 
quartz was displaced towards the curve for unirradiated silica. By x-ray, 
x-ray density, and optical methods it was later shown that extensive 
neutron irradiation reduced several phases of silica, including quartz and 
glass, to a common amorphous phase (Wittels and Sherrill 1954) and that 
this material reverted to «-quartz upon annealing at 930°c. 

The changes in density which accompanied the disordering process in 
quartz (Wittels and Sherrill 1954) followed a very peculiar behaviour, 
since, more often, the changes which are produced in solids as a result 
of high energy radiation follow a simple saturation curve when plotted 
as a function of irradiation dosage. In quartz, the relation of density 
change versus neutron flux exhibits a sigmoid curve with a maximum 
rate of density change occurring at ~ 6x 101° neutrons/cm?, and this 
rate approaches saturation with dosages in excess of 1-2 « 107° neutrons/cm?. 
Kinchin and Pease (1955) have suggested that the breakdown of the 
structure occurs when a sufficient concentration of interstitial ions are 
produced to render the structure unstable. The anisotropy of expansion 
(Wittels 1953, Johnson and Pease 1954) was thought by Wittels to be 
related to the open c-axis channels in the quartz lattice. Recently 
Klemens (1956) presented a theoretical treatment of this effect which 
was based on a mechanism of plastic flow that would result in an 
accelerated volume expansion rate. In a study of the optical effects 
in quartz resulting from neutron irradiation Mitchell and Paige (1956) 
interpreted their data as indicating point defects, namely oxygen inter- 
stitials and vacancies, and suggested that amorphous regions were formed 
in irradiated crystals as a result of local melting phenomenon. 

The purpose of the present paper is to discuss some recently obtained 
x-ray diffraction data on neutron irradiated quartz single crystals in 
relation to these hypotheses. 


§ 2. EXPERIMENTAL METHOD 


Single crystals of quartz were exposed to fast pile neutrons at 100° 
with dosage increments of approximately 1019 neutrons/em?. The 
integrated fast flux dosages are estimated to include all neutrons whose 
energies are in excess of 50 kev. Before irradiation and between 
successive periods of cumulative bombardment, the crystals were 
examined using X-ray diffraction and hydrostatic density methods. 

For purposes of making direct comparisons between the x-ray density 
data and the hydrostatic density measurements, cubes having an edge 


Structural Behaviour of Neutron Irradiated Quartz 1447 


length of 0-25 in. were cut from Y-cut bars so that {00-1} and {11-0} 
surfaces were coincident with cube faces. Precision lattice parameters. 
were measured by means of a spectrometer and geiger counter method 
employing a narrow incident beam (12 minutes total divergence). By 
the same technique, the shapes of several reflections were scanned by 
step counting in intervals of 0-01°(@). Hydrostatic density measurements 
were made on the same cube specimens so that a relation might be 
established between the calculated densities taken from the x-ray data 
and the densities measured hydrostatically. As will be shown presently, 
the x-ray diffraction measurements on these cubes are primarily 
concerned with the host crystal lattice and the effects which the 
irradiation-induced imperfections have upon that lattice. 

From the same crystal bars, thin plates varying between 0-003 in. 
and 0-007 in. were cut so that [00-1] and [11-0] directions were made 
perpendicular to the plate surfaces. These crystals were similarly 
irradiated and between successive stages of irradiation were examined 
by means of x-ray diffraction transmission methods that employed 
conventional single crystal film techniques on stationary and moving 
crystals at room temperature and at 78°K. When the appearance of 
“vitreous inclusions’ became apparent in film photographs, a geiger 
counter-step scanning method was used to determine the approximate 
shape of the diffuse x-ray scattering halo which is associated with the 
scattering of x-rays from an amorphous material. It will be seen that the 
x-ray diffraction transmission data furnish information concerning the 
irradiation-induced imperfections in addition to the long range order of 
the host crystal matrix. 

In a complementary study, polycrystalline samples of coesite, a high 
density crystalline silica (Coes 1953, Sosman 1954, Ramsdell 1955) were 
exposed to heavy dosages of fast neutron flux and examined with Debye— 
Scherrer x-ray diffraction methods. The purpose of these irradiations 
was to determine whether a possible structure dependence existed 
concerning the radiation stability of different phases of the silica system. 


§ 3. EXPERIMENTAL RESULTS 


3.1. X-Ray Diffraction and Density Measurements 


The results of the lattice parameter measurements are an extension 
of the earlier studies of the anisotropic expansion in irradiated quartz 
(Wittels 1953, Johnson and Pease 1954) and are shown in fig. 1. These 
measurements were taken from crystals irradiated with integrated fluxes 
only up to 8X 10!% neutrons/cm?, since beyond this dose, the reflections 
are too broad and diffuse for precision measurements. It is clear that 
the expansion in quartz resulting from neutron irradiation is anisotropic, 
and in addition, if one compares this irradiation-induced expansion 
(fig. 1) with the thermal expansion of unirradiated quartz (fig. 2) it is 
evident that a strong similarity exists between these two characteristics. 
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In both cases, the expansion perpendicular to the c-axis is significantly 
faster than the expansion parallel to the c-axis. Also, in both unirradiated 
and irradiated crystals, the rate of expansion Increases with either 
increasing temperature or increasing neutron dosage. 
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In the case of either unirradiated or irradiated quartz, it is apparent 
that the restoring force for atomic motion perpendicular to the c-axis 
is less than the same force parallel to the c-axis and this is verified by the 
anisotropy of expansion in both cases. This relationship could hardly be 
fortuitous and probably indicates that the irradiated quartz lattice is 
distended as a result of defect introduction in a manner which is consistent 
with the directional binding forces of the unirradiated crystal. 

Comparisons between the hydrostatic volume changes induced by 
irradiation and the volume changes calculated from the x-ray diffraction 
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data are shown in fig. 3. In this figure percentage changes in volume 
obtained with hydrostatic measurements are accurate to --0-02°/ while 
the same changes determined from x-ray data are accurate to -—-0-03°%. 
Within these experimental limits it is seen that for irradiation doses up 
to about 3x 10!* neutrons/em? the two kinds of measurements are in 
agreement, and further, that the changes in volume follow a trend that 
is similar to the thermal volume changes in unirradiated quartz (see 
inset fig. 3). The latter result is consistent with the observations of 
linear expansions previously given. The relationship between the unit 
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cell reciprocal volume as calculated from x-ray data, and the measured 
hydrostatic densities, provides evidence of a general nature concerning 
radiation damage in quartz. In fig. 4 it is noted that with irradiations 
up to ~ 3X 101% neutrons/em? the points fall on a straight line that 
is contiguous with the same straight line for unirradiated quartz, and 
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that at higher doses there is a marked deviation from this relationship. 
At low doses then, the molecular weights of irradiated crystals are 
identical to those of unirradiated crystals which contain three SiO, 
molecules per unit cell, and at higher doses, the increased concentration 
of defects having densities lower than those of the host crystal results 
in smaller calculated molecular weights. 
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It has been theoretically demonstrated (Miller and Russell 1952, 1953, 
Kshelby 1953) that the presence of Frenkel defects in an isotropic medium 
should result in an agreement between x-ray density and hydrostatic 
density values. Although this agreement is found in quartz irradiated 
with less than 3 x 101° neutrons/cm? is it not entirely clear that the theory 
is applicable to the highly anisotropic lattice of quartz. The steeper 
slope (fig. 3) of the volume change at high doses as determined by x-ray 
methods is in contradiction with Klemens’ model (1956) since ‘@ relief 
of shear strain through plastic flow of the crystalline regions would be 
expected to give rise to a saturation in the expansion of the host crystal 
matrix as determined from x-ray lattice parameters. 

Examination of the packing model of the quartz structure (Bragg and 
Gibbs 1925, Gibbs 1926) gives evidence suggesting that the configuration 
of interstitial space in the lattice may be related to the radiation damage 
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process in quartz. Nearly half of the total volume of the unit cell is 
unoccupied and more than two-thirds of this unoccupied space comprises 
the screw channels which are roughly parallel to the c-axis in a hexagonal 
array. It was suggested earlier (Wittels 1953) that these sites might 
serve as interstitial traps and thus qualitatively explain the anisotropy 
- of expansion. The oxygen ions which are situated immediately adjacent 
to the screw channels could be easily displaced to these sites by breaking 
only two bonds as against four bond ruptures for a silicon displacement. 
In addition, the large ionic diameter of oxygen would make it a more 
efficient defect for producing the anisotropic expansion. 

Except for line shifts, the x-ray diffraction patterns of crystals irradiated 
with less than 3 x 101° neutrons/cm? are relatively unaltered. With the 
limitations imposed by the total divergence of the incident x-ray beam 
(12 minutes), little or no line broadening was observed in this dosage 
range. This indicates that the expansion does not involve severe dis- 
tortion of the lattice at low doses, and that regions of highly disordered 
material of a vitreous nature are not present in concentrations sufficient 
for detection by this technique. It would be expected that the presence 
of large vitreous regions would set up inhomogeneous shear stresses 
(Klemens 1956) in the crystalline regions and result in the skewing of 
reflection peaks together with little or no peak shifting. For low doses 
at least, it is evident that this is not the case since the x-ray reflections 
remain sharp and are measurably shifted due to the lattice expansion. 
The x-ray diffraction patterns of crystals irradiated with neutron doses 
in excess of 6x 101° neutrons/em?, however, reveal considerable lattice 
strain as the Laue and Bragg reflections (fig. 5 (b), Pl. 59) are smeared 
and peak intensities are diminished. Klemens’ (1956) model predicts 
an extremely skew line broadening should result from inhomogeneous 
shear strain in the crystal matrix and experimental verification of this 
prediction is seen in x-ray diffraction lines of crystals irradiated with 
8x 1019 neutrons/em?. Figures 6 and 7 show that in the direction of 
greatest expansion (a-axis) the Bragg reflection is shifted and smeared 
to give a long tail of low intensity in the direction of increasing diffraction 
angle, while in the c-direction the reflection is shifted, but otherwise 
is much less affected. It is conceivable that the extreme skewing of 
(bh-0) reflections results from an inhomogeneous shear strain that is 
predominant in directions perpendicular to the c-axis. This idea is also 
consistent with crowding of interstitial ions into the open c-axis channels 
to cause the anisotropic expansion. In addition, the lack of comparable 
line skewing for (00-1) reflections supports this contention. The lattice 
parameters were calculated from the peak positions of the lines rather 
than the mean positions. Using the reflection peaks to calculate the 
x-ray volumes one finds that the ratio of x-ray volume to hydrostatic 
volume after irradiation with 8 x 1019 neutrons/cm? is 1-48. By employing 
the mean positions of the reflections the ratio becomes approximately 
1:14. This is in agreement with Klemens suggestion that line skewing 
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should account for a portion of the lattice expansion observed. It should 
be noted, however, that this extreme line skewing (fig. 6) is not observed 
in crystals irradiated with less than 7x 101° neutrons/em? and that 
lattice parameter measurements of these crystals are not exaggerated 
by the employment of reflection peaks. Although the calculated x-ray 
volume is reduced by employment of the mean positions of the reflections, 
it is evident that the x-ray volume is still greater than the hydrostatic 
volume after 8X 1019 neutrons/cm? and that there is no early saturation 
in the expansion of the host crystal as expected from Klemens’ model 
of plastic flow. 


Fig. 6 


BEFORE IRRADIATION 


oO 
E 
Psd 
=) 
> 
= 
oe AFTER 80 x 10'9 nem? 
= i /cm 
faa] ee 
a 
= 
> 
& 
w 
=a 
dij 
= 
2 
if 73 74 75 76 rae 78 79 80 


2 6 (deg) 


(22-0) of quartz single crystal irradiated with 8x 1019 neutrons/cm?, 


3.2. Diffuse X-Ray Scattering 
The studies of Berman (1951) suggest that disordered regions are 
produced in quartz by neutron irradiation and are observed in low 
temperature thermal conductivity measurements. In an optical investi- 
gation of irradiated quartz, Mitchell and Paige (1956) have suggested 
that amorphous regions are formed at all levels of irradiation Aeon h 
these defects are easily identified by X-ray measurements wii on es 
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irradiated with dosages in excess of 6X 10!* neutrons/em2, this does not 
preclude their presence in lower concentration at lower doses or in regions 
which still retain a semblance of order, or both. The changing nature 
of the disorder within these regions, and the complete vitrification to 
the amorphous state upon prolonged bombardment can be qualitatively 
observed by diffuse x-ray scattering measurements. 


Fig. 7 
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(00-3) of quartz single crystal irradiated with 8 x 101° neutrons/cm?, 


The diffuse x-ray scattering halo which is characteristic of amorphous 
materials becomes evident in single crystals of neutron irradiated quartz 
that still possess considerable long range order. As shown in figs. 8 (0), (c), 
Pl. 59, there is a broad diffuse halo with a peak position close to that for 
unirradiated vitreous SiO, that appears in the reciprocal lattice photograph 
of irradiated quartz single crystals. The shape and peak position of this 
diffuse scattering from irradiated single crystals is progressively altered 
(fig. 9) until the rate of density change becomes zero. Qualitatively, it 
is evident from these changes that the configuration of the atoms within 
the amorphous regions changes upon increased irradiation. The shift 
in the peak position is indicative of an alteration of the average inter- 
atomic distances within these disordered regions and is probably related 
to the constraints of the host crystal lattice as suggested by Klemens 
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(1956). Since the crystal matrix expands sharply during the development 
of these amorphous regions, it seems reasonable to assume that the 
increased openness of the host lattice facilitates the changing configuration 
of the disordered atoms and atom groups. It is also believed that ruptured 
Si-O bonds aid this rearrangement not only by furnishing a higher degree 
of freedom to [SiO,] tetrahedra, but also by providing the crystal lattice 
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with a driving force that tends to reduce the number of ruptured bonds 
by means of suitable rearrangements. The suggestion of a higher degree 
of freedom for [SiO,] tetrahedra has been proposed by Lukesh (1955) in 
studies on irradiated vitreous silica. In addition, the presence of 
ruptured bonds was also indicated by the magnetic susceptibility studies 
by Stevens (1955). 

Klemens (1956) has made estimates of the concentration of vitreous. 
material in quartz as a function of neutron flux by two separate methods. 


Structural Behaviour of Neutron Irradiated Quartz 1455 


These values were found to agree within a factor of 2. From the observa- 
tions of diffuse x-ray scattering in crystals studied here only qualitative 
estimates can be given concerning the fraction of vitreous material 
present in irradiated single crystals. These estimates tend to agree with 
the lower values presented by Klemens and indicate that 10-15%, of the 
material is vitrified after irradiation with 5x 1019 neutrons/cm?2. 

The Laue photographs of crystals irradiated with 7 x 101° neutrons/em? 
(fig. 5 (6), Pl. 59) exhibit a diffuse scattering pattern of spots and streaks 
in addition to the halo that is associated with the amorphous regions 
of the crystal. This complex diffuse scattering pattern becomes even 
more evident when monochromatic radiation is used since Laue and Bragg 
diffraction is then suppressed. With single crystals oriented so that the 
incident monochromatic beam is parallel to the c-axis, a pattern of 
scattered radiation results which is outside the directions of selective 
reflection and which exhibits a hexagonal symmetry that is related to 
the pseudohexagonal symmetry of the quartz lattice. 

In fig. 10 (a), Pl. 59, the six intense reflections of the outer hexagonal 
shell and the six strong inner spots are Bragg reflections while the 
remaining spots and streaks occur outside the zones of selective 
reflection. The appearance of the Bragg reflections is primarily due to 
the irradiation-induced distortion of those particular planes which happen 
to lie very close to the proper position for satisfying the Bragg equation. 
To a very minor extent this selective reflection is also due to a small 
divergence of the incident monochromatic beam. 

The diffuse spots and streaks that fall outside the directions of selective 
reflection are remarkably similar to the anomalous scattering of a thermal 
origin reported by several workers (see Lonsdale 1942-43). It is unlikely 
that a thermal origin can be employed to describe the formation of this 
diffuse scattering in neutron irradiated quartz, however, and it is more 
probable that it is predominantly due to the displacement of the atoms 
from their mean positions by the presence of defects in the lattice. In 
this regard it should be emphasized that similar diffuse patterns, as in 
fig. 10 (a), Pl. 59, are conspicuously absent for unirradiated crystals 
under the same conditions of exposure time and monochromatic radiation. 
Under these conditions the photographs do not contain the Bragg spots 
of fig. 10 (a), Pl. 59, because of the lack of distortion of the unirradiated 
lattice, and the diffuse pattern is not even weakly evident. In addition, 
since the intensity of diffuse scattering from the irradiated crystals is 
hardly affected by cooling the crystals to 78°K, this is further evidence for 
ruling out a thermal origin in quartz which has a characteristic tempera- 
ture, 9, of approximately 300°K based on acoustical mode frequencies. 
It appears from the oscillator data of Johnson and Pease (1954) that 
changes of less than 1° are produced in the elastic constants, and it is 
therefore inconceivable that these small changes could alter significantly 
the thermal diffuse scattering. From these results, it appears that the 
anomalous diffuse x-ray scattering observed in heavily irradiated quartz 
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is not of thermal origin but probably results from the inhomogeneous 
strains induced by the bombardment-produced defects. 


3.3. Neutron Irradiation of Coesite 


The classical studies of Bragg (1930) demonstrated the unique role of 
the [SiO,] tetrahedral group in silicate crystals. In the crystalline phases 
of silica these structural entities are joined together by a sharing of all 
tetrahedral corners to form a rigidly bonded network which simply means 
that each oxygen ion forms a bridge between adjoining tetrahedra. The 
same structural unit is an integral part of vitreous silica as well (Warren 
et al. 1936), but in this case the tetrahedra are irregularly linked to form a 
random network configuration. From a structural viewpoint it is there- 
fore evident that the basic difference between silica solids is primarily due 
to the Si-O-Si bond angle relationships and bond distances which in large 
part determine the closeness or openness of atom packing. 

In quartz, low tridymite, low crystobalite and vitreous silica, the atom 
packing is markedly open and results in 50°, or more void space in all 
of these solids. Since these materials are reduced to a common amorphous 
solid following prolonged neutron bombardment (Wittels and Sherrill 
1954) it was of interest to determine whether the radiation damage 
process in silica structures was related to the openness of structure. For 
this purpose, polycrystalline specimens of coesite, a high density synthetic 
silica, were exposed to 2 x 10° neutrons/cm?, the dosage which is sufficient 
to vitrify three ‘open’ silica structures. Before and after irradiation 
Debye-Scherrer x-ray diffraction patterns were taken to ascertain whether 
any structural changes had occurred in the material as a result of the 
neutron bombardment, and it was determined that no appreciable damage 
had been frozen-in since diffraction patterns remained relatively unaltered. 
Lattice expansions in excess of 0-1°%, were not observed and the diffraction 
lines were not noticeably broadened. 

Although the structure of coesite (Ramsdell 1955) is not clearly under- 
stood it is reasonable to assume that the crystal is composed of linked 
[SiO,] tetrahedra similar to the natural silica forms. The essential 
difference between coesite and these natural forms of silica is therefore 
displayed by the packing of atoms in the structure. The synthetic 
material has a density of 3-01 g/cm? and must be considered close-packed 
when compared to the other phases of silica. Presumably the binding 
energies in coesite are not considerably greater than the binding energies 
in the natural silica phases, and if this is true, the production of displace- 
ments by knock-on collision should not be significantly different for either 
case. And since the basic physical differences between these structures 
concerns the packing of atoms it seems reasonable to believe that the 
radiation damage process in silica structures is primarily concerned with 
the ease or difficulty with which [SiO,] tetrahedra may be reoriented. 
In the open silica lattices it is conceivable that these normal voids furnish 
ideal trapping sites for interstitial ions while the acompanying rupturing 
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of bonds assists in ‘loosening’ the lattice to make the structure more 
susceptible to further damage. This process is necessarily obscure but 
it is suggested that the radiation stability of coesite offers evidence 
supporting the suggestion that both the rate of damage and the location 
of lattice defects have a structure dependence in neutron irradiated 
silica solids. 


§ 4. Discussion oF RESULTS 


The structural effects of the neutron irradiation of quartz with doses. 
up to 1-2 x 10° neutrons/em? follow at least two processes : (i) an aniso- 
tropic expansion that is essentially elastic at doses less than ~3 x 101° 
neutrons/cm*, and (ii) a more complex anisotropic expansion at higher 
doses. Associated with the first process is a direct relation between 
x-ray densities and hydrostatic densities. With the second process one 
observes: (a) inhomogeneous shear strain in the crystalline matrix 
perpendicular to the c-axis, (b) a growth and changing configuration of 
amorphous regions within the host lattice, and (c) a diffuse x-ray scattering 
pattern indicative of inhomogeneous lattice strains that result from the 
presence of defects. Coesite, a high density crystalline silica is stable 
under neutron irradiation doses which completely disorder quartz. 


4.1. Effects after Doses Less than 3 x 1019 


A model of the radiation damage process in quartz is developed on the 
basis of point defects and slightly disordered regions in crystals irradiated 
with ~101%3 x 101° neutrons/cm?, in the range where the lattice dis- 
tentions appear to be elastic. This analysis appears to fit the annealing 
studies of Berman (1951), Johnson and Pease (1954) and Mitchell and 
Paige (1956) on crystals irradiated with ~1 101% neutrons/cm?.  Ber- 
man found that the complex defects annealed more rapidly than the 
simple defects and this is consistent with a picture of only slightly dis- 
ordered regions plus simple defects. Primak (1955) argues that the 
absence of any stored energy release is proof for the non-existence of inter- 
stitial-vacancy defects in neutron irradiated quartz. It is conceivable, 
however, that a possible stored energy release would remain undetected, 
either because of its small magnitude, or by virtue of the fact that a slow 
release over a broad spectrum of temperature would be difficult to deter- 
mine experimentally. The presence of point defects, first observed by 
Berman, tends to be confirmed by the annealing recovery of the lattice 
parameter and piezoelectric properties (Johnson and Pease 1954) as 
well as the optical properties (Mitchell and Paige 1956) in irradiated 
crystals. More recently, magnetic studies (Weeks et al. 1957) have also 
revealed the presence of point defects in irradiated quartz. It is suggested 
that the essentially complete recovery of order by thermal annealing 
would be unlikely if all the defects were of large and complex configuration 
since the rigidity and short range character of binding in the silicon— 
oxygen system is not amenable to an annealing process that requires. 
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mobility. In the present study the lack of distortion seen in x-ray diffrac- 
tion patterns of crystals irradiated with less than 3x 101° neutrons/em? 
also supports these conclusions. ne 
The displacement calculations of Mitchell and Paige (1956) indicate 
that on the order of twice as many oxygen ions would be displaced by 
knock-on collision than silicon ions. Although this may not be con- 
sidered a large enough factor to insure the suggestion that oxygen ions 
form the predominant interstitials, it is still evident that the large ionic 
diameter of oxygen makes it a very efficient agent for inducing expansions 
of the lattice. Tentatively, therefore, the present author is in agreement 
with the suggestion of Mitchell and Paige concerning the presence of 
oxygen ion interstitials and vacancies in neutron irradiated quartz. 


4.2. Effects after Doses in Excess of 3 101° neutrons/em? 


The effects of irradiation with doses in excess of 3 x 101° neutrons/cm? 
are exceedingly complex and still not well understood. Nevertheless, 
the experimental results at these high doses are compatible with the 
conclusions given concerning crystals having suffered less damage. The 
extreme skewing of the (22-0) reflection compared to the rather slight 
broadening of the (00-3) reflection indicates that an inhomogeneous shear 
strain is present in crystals irradiated with 8 x 101% neutrons/em?. In 
part, this line skewing is in accord with the prediction of Klemens (1956) 
but it is even more consistent with the anisotropic expansion of the lattice 
and is propably due to the crowding of interstitial oxygen ions into the 
open c-axis channels. 

It is inherent that ruptured silicon—oxygen bonds are associated with 
the production of point defects and ‘ vitreous’ regions since the un- 
irradiated perfect crystal is composed of a rigidly bonded network of 
[SiO,] tetrahedra. It is difficult to conceive of the disordering of four 
phases of silica (Wittels and Sherrill 1954) to a common amorphous phase 
without the action of some rearrangement process that would promote 
the formation of a common silicon-oxygen configuration. This process 
is probably assisted by the presence of ruptured silicon—-oxygen bonds. 
It is also suggested that the thermal recrystallization of quartz from 
heavily irradiated silica phases (Wittels and Sherrill 1954) also is aided 
by the presence of broken silicon-oxygen bonds since, in either case, the 
processes would be promoted by the presence of ruptured bonds which 
would result in a loosely held array rather than a rigid network. 

The appearance of a diffuse x-ray scattering halo in diffraction patterns 
of heavily irradiated quartz (7 x 1019 neutrons/cm?) confirms the presence 
of amorphous regions in crystals that still possess considerable long range 
order. With increasing dosages of irradiation the growth and modifica- 
tion of these amorphous regions is indicated by the change in peak 
position, shape, and relative intensity of the diffuse halo. The exact 
nature of this phenomenon is not clear, but it is evident that the average 
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interatomic distances and atomic configuration within these regions 
changes as the host lattice rapidly expands. Qualitatively, the shift 
of the diffuse peak to higher values of sin 6/A means that the average 
interatomic distances within the amorphous regions is smaller than that 
found in unirradiated vitreous silica. Perhaps it is unusual to find that 
the specific volume of these amorphous regions does not increase sharply 
as the host lattice rapidly expands but since the density of irradiated 
fused silica increases upon neutron irradiation (Primak et al. 1953) it 
appears that the process is not one of simple expansion. It is probably 
also related to the suggestion by Klemens (1956) that the volume ex- 
pansion of these disordered regions is constricted by the surrounding 
crystalline matrix. The diffuse scattering studies of Simon (1956) on totally 
disordered quartz indicated that the amorphous material showed a shorter 
Si nearest Si neighbour distance, and a smaller Si-O-Si bond angle than 
that given in the model of silica by Warren et al. (1936). In addition, his 
study showed a wider distribution of both O-O distances and more 
distant Si-O coordination distances which might be considered as evidence 
for randomly placed point defects and ruptured Si-O bonds. 

There is no positive evidence in the present investigation which would 
indicate that a thermal spike mechanism (Primak 1955) is the cause for 
the irradiation-induced disordering of quartz. Moreover, it is difficult 
to see how the defects produced by Seitz’ knock-on mechanism could 
be distinguished from the disordering due to thermal spikes (Seitz 1949, 
Brooks 1950) or displacement spikes (Brinkman 1954). If an oxygen 
ion were displaced to an interstitial position by the knock-on process, 
in the vicinity of the oxygen ion vacancy the first oxygen coordination 
shell for two adjoining [SiO,] tetrahedra would be disrupted and this 
would probably be followed by a rearrangement of these shells to better 
satisfy the binding requirements. The coordination shell rearrangements 
would result in an alteration of the Si-O-Si bond angles as well as bond 
distances in a small region about the vacancy and would cause a sharp 
decrease, locally, of short range order, and a much smaller decrease in 
long range order. In addition, a distortion would be affected in the 
region of the interstitial oxygen ion which would also give rise to local 
disorder. In a necessarily naive sense, these disordering effects in a 
refractory crystal like quartz would be much like a ‘ vitrification, > in 
local areas the sizes of which would depend upon the local concentration 
of interstitials and vacancies. Where displacements are more closely 
spaced near the end of the track of the primary recoil ion these * vitrified ’ 
regions would occupy a larger volume, as suggested by Klemens. 


§ 5. CONCLUSIONS 
X-ray diffraction studies of neutron irradiated quartz indicate that at 
least two damaging processes occur with doses up to 1-2 x 10° neutrons/ 
em?, Crystals irradiated with less than ~3 x 101° neutrons/cm? exhibit 
an anisotropic lattice expansion which is similar to the thermal expansion 
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of unirradiated quartz, and it is suggested that this is related to the 
directional binding forces of the unirradiated crystals. In this dosage 
range, the volume increase is the same, as determined by x-ray and 
hydrostatic methods, which might be interpreted as indicating the 
presence of point defects. Except for line shifts, x-ray diffraction 
patterns from single crystals show no observable distortions and indicate 
that the expansion is essentially elastic. Although it is felt that point 
defects and ‘vitreous’ regions are produced by knock-on collisions 
at all doses, it is believed that the total disorder of the complex defects is 
relatively slight in crystals irradiated with less than 3 x 101° neutrons/cm?. - 

At doses in excess of 3x 101° neutrons/cm? the damaging process is 
more complex since the volume change indicated by X-ray measurements 
of the host lattice becomes larger than the bulk volume change measured 
hydrostatically. The inhomogeneous shear strain induced by the large 
anisotropic expansions after 8 x 1019 neutrons/cm? is seen in the skewing 
of the (22-0) reflection and reveals the anisotropy of the damaging process 
after high doses. Amorphous regions enmeshed within the host lattice 
are identified in crystals that still possess long range order after irradiation 
with 7x 101° neutrons/em?. The diffuse x-ray scattering halo associated 
with this amorphous material changes peak position, shape, and in- 
tensity as the host lattice is distended and strained, and is probably 
due to the rearrangement of the atoms within the disordered regions 
as the host matrix expands. In addition, a second kind of diffuse scatter- 
ing appears which cannot be accounted for on the basis of a thermal 
origin but is evidently related to scattering as a result of inhomogeneous 
strains in the damaged lattice. 

The stability of polycrystalline samples of coesite, a high density 
synthetic silica after irradiation dosages sufficient to completely disorder 
the natural ‘ open ’ phases of silica suggests that the rate of damage has a 
structure dependence in neutron irradiated silica solids. In neutron 
irradiated quartz, there is some evidence that the lattice defects may be 
structure dependent as well, but this is less certain for the other phases 
of silica and coesite. 
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§ 1. INTRODUCTION 


Crystats which conduct electricity by movement of the sons in the 
crystal have proved rather favourable subjects for the study of lattice 
defects such as interstitial ions, vacancies and foreign ions. In particular, 
much is now known about their concentrations, their mobilities and their 
interactions with one another. With the aim of extending this knowledge, 
particularly of the kinetic properties of the lattice defects, we have studied 
the thermoelectric power of ionic conducting crystals. Previous work on 
this topic has been almost entirely from the point of view of the theory 
of irreversible thermodynamics (de Groot 1951, Holtan 1953). We have 
approached the subject theoretically by making use of the models of 
Frenkel and Schottky lattice disorder and have obtained results which 
are consistent with the thermodynamic analysis but which allow more 
detailed predictions to be made. 

We shall here confine ourselves to a crystal MX fitted with electrodes 
of metal M, the experimental arrangement being as shown schematically 
in fig. 1. The thermoelectric potential is the difference in potential 
V,—V,4, made of up five parts: (1) the ‘homogeneous’ potential 
difference Vp,p—V_, due to the temperature difference between P and A ; 
this term is small by comparison with the others and will be dropped ; 
(2) the contact potential difference Vo—Vp= —O@(T); (3) the homo- 
geneous potential difference V,—V, due to the temperature gradient 
in the salt; (4) the contact potential difference V,—Vp=@(T'-+ JT) ; 
(5) the homogeneous term in the electrode and wire which we shall again 
neglect. The total is then 


Vu—Va=(Va—Vq) + G(T + AT)— G(T), 

The quantity (V,—V,)/AT is the thermoelectric ‘power’ 6 and is 
composed of a homogeneous part Oomn=(Vr— V.)/AT and a hetero- 
geneous part @ne-=46/AT. We shall discuss these two components 
separately. In this discussion it will be convenient to have a particular 
system in mind; we shall therefore suppose we are dealing with a 
cationic conductor displaying Frenkel disorder, i.e. containing mobile 

ee ee 
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interstitial cations and cation vacancies. This is typified by AgBr or 
AgCl. In addition we shall deal with crystals containing foreign ions 
of a different valency, e.g. Cd?+ in AgBr. The importance of this case 
comes from the creation of one Ag+ vacancy for every molecule of CdBr, 
incorporated in the AgBr crystal; the vacancy concentration can 
therefore be altered at will. 


Fig. 1 


Ul et A 


A B 


Schematic diagram of the experimental arrangement for the measurement 
of the thermoelectric power of an ionic conductor MX fitted with 
electrodes of metal M. In the equations of §2, the currents 7 and the 
spatial coordinate x are referred to left to right as positive direction. 


§ 2. HomocEnEous Power, @nom 


The thermoelectric potential is measured with no current flowing, so 
we therefore write down phenomenological equations for the current 
densities of the two charge carriers (interstitials and vacancies) 


co adtl 
j=L,| —(erad pe)y—e grad V4 |: elie 


1 9P : 
ju=L,| —(grad 1,)y be rad y- qs |, er) 
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and determine the power ?rom=(grad V)/(grad 7’), by: setting the total 
electric current equal to zero, Le. 
Ji —J y=. utes dol III 2a) 


The coefficients L, and L,, are proportional to the mobilities of interstitials 
and vacancies respectively and q,* and q,* are ‘ heats of transport’, — 
i.e. coefficients of proportionality between the defect current densities 
and the ‘ reduced ’ heat flow under isothermal conditionst. The chemical 
potentials », and yw, are defined as the change in the Gibbs free energy 
of the crystal on adding respectively one interstitial and one vacancy. 
From the statistical theory of Frenkel disorder we obtain the expressions 


pmgb kt In (5), . ag eee 
via 

= = Fae « 5: De 

cs +k In(“2), aa (2.5) 


where g, is the work required to bring a cation from a state of rest at 
infinity into a particular (but arbitrary) interstitial position in the 
crystal (constant 7’, P). Similarly g, is the work required to remove a 
normal ion to a state of rest at infinity. N and N’ are the numbers of 
normal and interstitial sites per unit volume and n,; and n, are the 
numbers of interstitial ions and vacancies per unit volume. Equations (2.4) 
and (2.5) allow us to write eqns. (2.1) and (2.2) in terms of grad n; and 
grad n, and these quantities are related to each other and to grad 7 
through the condition of local thermodynamic equilibrium which gives 
the familiar equation 


1 
NN, »=NN’ exp | ai (a9) | : , Be one 
This procedure gives for Onom the expression 


grad V __ kf gradn, , ((9,.*+h)é—Ga*/€} 
grad T en, grad T' eT (€+4/€} ; 


Shom = (2.7) 


where 


is the heat of formation of a Frenkel pair, and where ¢ is the ratio of 
interstitial mobility to vacancy mobility and € the vacancy concentration 
relative to its value in the pure crystal at that temperature (this may 
differ from unity owing to the presence of foreign ions of different valency) 


+ Holtan. (1953). The movement of a defect from one site to another is a 
thermally activated process and requires that energy be given to the defect 
at the initial position and given up again at the final position. There is 
therefore a heat flow associated with a flow of defects, additional to that 
resulting from the dependence of free energy on concentration (which in (2.1) 
and (2.2) gives the grad » terms). It is this additional part which corresponds 
to the ‘ reduced’ heat flow of the thermodynamic treatment. ~ 
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§ 3. INHOMOGENEOUS POWER, Ohet 


This is found from the contact potential 6(7') by differentiation with 
respect to 7’. The contact potential is determined from the condition 
that the electrochemical potential of M+ be the same in the salt and in 
the metal electrode, i.e. 


l 
0(T)= _ {u;,—u(M+t in M)}. ar ey eek 
Therefore 
pee ain 3.2 
het= 7 am + 78 Meal Jag we ee at Ree aa (ose) 


where s(M~* in M), the partial entropy of M~* ions in the metal, can be 
identified with the entropy of the metal since the free electrons make 
only a small contribution. The derivative du,/0T is found from (2.4) ; 
we get finally 


es ce Hse kT gradn, | 1 as k 
Cpe i + 7m (3) 4 en, orad i t 7 (Mt in M)— 7 a (3.3) 
where s;=—0g;,/0T, and nj is the value of n; for the pure crystal. 


§ 4. ToraL Power, @ 
This is obtained by adding together (2.7) and (3.3). We obtain 


— {e* +hé—4b9,*/€} & 1 Aran ret (5) 
0(é)= E+ E/E} ving - 7 S(AE* in M) + ~In y’)* 


(4.1) 


We shall not discuss here the value 6(1) of the power for the pure crystal 
(Howard and Lidiard 1957) but will consider the difference in power 
between a crystal containing a certain concentration of foreign ions 
(giving € different from unity) and a pure crystal. From (4.1) it follows 
that ae 
(?—1) bai" +a." +h) & ry 
~ (2446) @+DeP mae es) 


The relative vacancy concentration € may be obtained from a theoretical 
discussion of the influence of foreign ions on the defect concentrations or, 
in favourable cases, can be derived directly from experimental results 
on ionic conductivity, e.g. in AgBr and AgCl (Teltow 1949, Ebert and 
Teltow 1955). In these systems the conductivity isotherms are described 
by the equation 


Ad= 6(€)—A(1) 


hee, Sees | (4.3) 
99 1+¢ 
where o is the conductivity of the impure crystal, o9 is the conductivity 
of the pure crystal and ¢ as before is the ratio of mobilities. For both 
AgBr and AgCl the curves of c/o) against impurity concentration show 
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minima, the positions of which enable ¢ to be determined directly 
(Teltow 1949). It is then possible to use the experimental conductivity 
isotherms to give € as a function of impurity concentration c and then 
to calculate 6(€)—0(1) as a function of ¢ using (4.3). Typical results 
for AgBr containing CdBr, are shown in fig. 2. There is one unknown 


Fig. 2 
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° 
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ee] 
a 
Ww 
= 
9 eee 
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MOLAR PERCENTAGE OF Cd Br, 


The een power of AgBr containing CdBr, relative to the power of 
ae gf rat 175°c. These curves have been constructed by assuming 
va Bs for the parameter (q;*+q,*-+h) and by finding € from the 
ee eee pe ae of Teltow (1949). It may be noted that 

e power 0 75°C is fe i 
en pure AgBr at 175°c is —1365 pv/°o (Patrick and Lawson 


parameter (q¢,*+¢,*+h), whose value should be given rather accuratel 
by measurements of thermoelectric power as function of concentrati oy 
For reasonable values of this parameter the changes ares 


cad in power are lar 
and have a very characteristic variation with c oy 
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More generally the dependence of 6()—6(1) on ¢ can easily be discussed 
qualitatively since € is a monotonically increasing function of c. We find 
that whenever (q7;*+¢,*+h)/kT >(1+¢)?/26 then 6(£€)—0(1) initially 
Increases with c before passing through a flat maximum as in fig. 2. 
When (9,*-+-¢,*+h)/kT <(1+¢)?/26 then 6(€)—6O(1) initially decreases 
with increasing c; a minimum and a maximum occur at €— and &, 
given by 

4° =(e—NAd{(e—1P—1})™, 


where «=(q;*+9,*+A)/kT. &, and € are real if «>2. For «<2, 
6(€)—6(1) is a monotonically decreasing function. In assessing the 
magnitude of (¢;*+-q,*-+-h) it may be noted that q¢,* is a negative quantity 
whose magnitude is probably close to the activation energy for vacancy 
movement ; q¢;* is positive and probably rather smaller than the activation 
energy for interstitial movement. 


§ 5. CONCLUSION 


In conclusion we propose that measurements of thermopower of 
‘doped ’ ionic conductors will (a) provide a rather general and unambiguous 
check on the theory of lattice disorder in these crystals, and (b) provide 
a good value for q;*+q,* since can be determined independently. 
The values of ¢;* and q,,* are closely related to the details of the movement 
of lattice defects from one position to another. At present there exists 
only a very qualitative theory of q,* and q¢,*, and the development of 
crystal dynamical theory to improve on this situation would seem very 
desirable. Measurements of thermoelectric power in the systems we 
have described might provide a good stimulus for that development. 
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Lifetime of the 3-56 MeV State of Lithium 67; 
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RecentLy Browne and Bockelman (1957) obtained the first precise 
location of the second excited state of ®Li—3-560+.0-006 Mev—from 
an average of several observations on the reactions SLi(p, p’)®Li and 
*Be(p, «)*Li. They compared this value to the previously determined 
(Mackin 1954) gamma-ray energy of 3:572--0-012 Mev and concluded 
that the agreement in energy indicates that the 3-56 Mey level is relatively 
long-lived (comparable to the stopping time of the ®Li nuclei)s. The 
3-56 Mev level of ®Li is almost certainly the analogue of the *He ground 
state and thus has J7=0+, 7=1. In L-S coupling the M1 lifetime of the 
3-56 Mev level of *Li is calculated to be 5-5 10-17 sec. The lifetime can 
also be calculated using the M1 lifetime formula given by Sherr et al. (1955) 
which is valid for p? and p-? configurations in any degree of intermediate 
coupling. [ff the *He beta decay ft value and the ®Li magnetic and 
quadrupole moments are used to estimate the wave functions of the ®Li 
levels involved, it is found that the lifetime calculated from the formula 
of Sherr et al. differs negligibly from that obtained in pure L-S coupling. 
Therefore, a Doppler shift less than the maximum possible would indicate 
either that the transition were not M1 and thus the 3-56 Mev level were 
not 0+ as expected, or that the M1 lifetime was in serious disagreement 
with theoretical predictions. It should be noted that the above quoted 
energy measurements do not demand a relatively long-lived level since 
the quoted errors in the two measurements are also consistent with the 
full Doppler shift of 26 kev. However, in view of the serious consequences 


{ Communicated by the Authors. Work performed under the auspices of the 
U.S. Atomic Energy Commission. 
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§ The gamma-ray energy was obtained using the ®Be(p, «)®Li reaction at the 
2:56 Mev resonance. Stoltzfus et al. (1956) have shown that the 3-56 Mev 
gamma-rays following this reaction have an isotropic distribution and that the 
alpha—gamma correlation in the plane containing the beam is isotropic to 4%. 
Since the reaction proceeds through the 1°B 8-89 mev level, the angular 
distribution of the excited Li nuclei must be symmetric about 90° in the 
centre-of-mass system, and thus the full Doppler shift (defined as the shift 
expected for a lifetime short compared to the stopping time of the excited nuclei) 
is given by 4H=E{(v/c) cos 0}, where v is the velocity of the centre-of-mass 
system. The full Doppler shift expected for Mackin’s measurement (made at 0° 
to the proton beam) is 26 kev. 
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if, in actual fact, the 3-56 Mev level were found to be relatively long-lived, 
it was felt worthwhile to investigate the lifetime by a Doppler. shift 
measurement using the °Be(p, «)®Li reaction. 

Three different targets were used in this experiment. A thick beryllium 
target, a thick copper—beryllium target containing 2°%, beryllium, and a 
460 A beryllium target evaporated onto a 0-00005 in. nickel foil. For the 
thick target measurements, protons of initial energy 2-67 Mev were used. 
This energy was picked so that an integration of the measured cross 
section (Day and Walker 1952) along the proton path gave an effective 
proton energy corresponding to the peak of the 2-56 Mev resonance. For 
the thin target measurements, the beryllium was bombarded through the 
nickel foil and the proton energy was adjusted for maximum gamma-ray 
intensity. The gamma-rays were detected using a 2 in. right cylinder of 
Nal(Tl) mounted on a DuMont type 6292 photomultiplier which was 


Paes NT! vi. alae Set clita ee oe 
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CHANNEL NUMBER 


The full-energy-loss peak of the Li 3-56 Mev gamma-ray produced by the 
*Be(p, «)®Li reaction. The spectra were recorded by a 2 in. x2 in. 
Nal(TI) crystal and a 100 channel analyser with a bias of 400 channels 
at angles of 0° and 160° to the proton beam. 


picked for its non-rate-dependent gain. The crystal was irradiated with 
its centre 14-7 cm from the target, at 0° to the proton beam and at back- 
ward angle of 160°. The gamma-ray spectra were recorded using a 
“hundred channel analyser with a bias of 400 channels. Only the full- 
energy-loss peak was investigated. A typical set of spectra obtained from 
the thick beryllium target is shown in the figure. 
To obtain the pulse height of the full-energy-loss peak, the background 
was subtracted and the centre of gravity of the remaining peak was found 
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numerically. This method was possible because of the large number of 
experimental points in each peak. The gain of the system was monitored 
and corrections were made for shifting of the gain with time and position 
by comparison with the gamma-rays from the decay of ®°Co in a manner 
similar to that used by Wilkinson (1957) in his measurement of the 
Doppler shift of the first excited state of 1B. In all, 16 spectra were 
obtained in one sequence from the Be target, 44 spectra were obtained in 
four sequences from the Cu-Be target, and 27 spectra were obtained 
in two sequences from the Ni-Be target. The mean absolute shift in the 
pulse height of the full-energy-loss peak obtained from these three targets 
is shown in the second column of the table. 


Results of the Doppler Shift Measurement on the 
3:56 Mev Gamma-Ray from ®Li 


Measured absolute : | 
Target shift (in channels) Percent shift 
Be 5-83 10:2 1-43 +0-05 
Cu-Be 5-54-0-1 1-40-+0-03 
Ni-Be 5-62-40-1 1-42 +0-03 


The uncertainty given is the standard deviation. Because the 
amplifier-analyser system was found to be slightly non-linear, the 
analyser region between 400 and 500 channels was calibrated from the 
60Co 1:17 and 1:33 Mev and the 72Na 1-27 Mev gamma-ray lines. In 
order that the pulse heights of these gamma-rays would be comparable to 
the beryllium gamma-ray, the beryllium spectra were taken using a pre- 
cision attenuator which was removed when the radioactive spectra were 
taken. The attenuator was calibrated by comparing pulse heights of 
different gamma-ray lines, obtained with and without the attenuator, in 
a region of pulse height where the system was linear. The Doppler shifts 
expressed in percent, which are given in the last column of the table, were 
then obtained by converting the channel shifts to Mev by means of the 
calibration and dividing by 3-56 Mey. The errors here given include the 
uncertainty in the calibration as well as the uncertainties in the channel 
shifts. The observed shifts are to be compared to the calculated full 
shift of 1-44°%, (which includes a 0-5°% correction for finite counter size). 

It is seen from the table that the measured Doppler shift is constant 
within the experimental error over a wide range of stopping material 
density, and thus of *Li stopping time, and agrees, within the experimental 
error, with the calculated full Doppler shift. Therefore it is concluded 
that the lifetime of the 3-56 Mev level is short compared to the stopping 
time of the ®Li nuclei. The sharpest upper limit for the lifetime of the 
3-56 Mev level can be obtained from the Cu-Be measurement since the 
stopping time of the ®Li nuclei is shortest in the copper target. The 
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stopping time of the *Li nuclei in copper for the pertinent velocity region 
was extrapolated from data on the slowing down of 7Li nuclei in copper 
measured for slightly higher velocities (Devons and Towle 1956), and 
from the slowing down of 7Li in gases (Allison and Littlejohn 1956) 
measured in the pertinent velocity region. Before constructing a velocity 
versus stopping power curve covering the pertinent velocity region, the 
gas stopping power data was lowered by approximately 20%. This 
correction for the difference in the stopping powers of metals and gases 
was estimated from the published (Allison and Warshaw 1953) dependence 
of the stopping power of gases and solids for alpha particles which covers 
a wide range of velocities and densities. From these data it is found that 
a Doppler shift of 1:30°, for the copper target would give a limit 
T 42<3X 10-14 sec. Since 1-30°% is three standard deviations below the 
measured shift of 1-40--0-03°%, it is concluded that the lifetime of the 
SLi 3-56 Mev level is less than 3x 10-14 sec. This limit is sharp enough to 
give strong evidence against an E2 transition. 

The Weisskopf formula would present 9x 10~! sec if the transition 
were E2. The worst possible (and hardly imaginable) inhancement of 
this transition due to collective motion would arise if all three protons 
contributed to the transition, in which case the E2 lifetime would be 
10-13 sec. Therefore, it is concluded that the transition is predominantly 
M1 in character in agreement with the results (Mackin 1954) obtained from 
the internal pair spectrum of the *Li 3-56 Mev gamma-ray. 
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Scattering of Polarized Positrons by Polarized Electrons 
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ABSTRACT 


The cross sections for the scattering of polarized positrons and electrons 
by polarized electrons have been calculated. 


$1. INTRODUCTION 


THE usual methods of measuring the spin orientation in beta rays depend 
on scattering by heavy nuclei. The spin effect in this phenomenon, how- 
ever, is quite small in the case of positrons (Massey 1943). As an alterna- 
tive method for detecting the spin-orientation, Dicke (1950), Frauenfelder 
(1957) and others have suggested methods depending on scattering by 
polarized electrons. 

Since there has been some doubt as to whether the discrimination would 
be big enough, especially with positrons, we have calculated the cross 
sections for each kind of polarization and for bothet ande~. The electron— 
electron cross section given by Moller (1932) and the positron—electron 
cross section given by Bhabha (1936) refer of course to unpolarized 
electrons; their results follow from ours after averaging over the spins. 
The calculation is a straightforward one (Jauch and Rohrlich 1955) ; it was 
in fact carried out using eigenfunctions of op in the standard representation 
rather than by the usual method of traces. 

One or two errors were found in the literature; in particular we are not 
aware of any reference in which the last two terms of Bhabha’s formula 
are quoted correctly. 

§2. RESULTS 

A positron with total energy ymc? collides with an electron at rest, whose 
spin in the direction of motion of the positron is +}. The spin of the 
positron in this direction may be (@) +4 or (b) —}. The positron loses a 
fraction w of its kinetic energy (y—1) mc?, and is deflected through an 
angle \ in the laboratory frame of reference ; @ is the scattering angle in 
the centre-of-mass system of coordinates. Thus 


w==4. (108.0). 924 
and 


tan 40=[4(y-+1)]* tan. et) 


{+ Communicated by the Author. 
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The differential cross section is then 


re? . nit GES 
1S yg We eg a (OR RS Wc 
Sy)" Hynes [1+3(y—1) sin? Tp ete sinAdd (3) 


where r')”=(¢?/mc*)? and.Q is given by the following formulae : 


(4a) 


ed Sypll py ts 
O) Oe Eat Te TFL aa et doth | 
(45) 


The average of these is equivalent to Bhabha’s formula (Bhabha 1936, 
Ashkin et al. 1954); it is 
Set eh oa eee ci near dy is ws Lue BG 
ET en 2 Ve 
If the incident particle is a negative electron, the values of Q to be inserted 
in (3) are: 


(5) 


2  (y—1)?+2y(3y+1) cos? 6+ (y?—1) cos! 6 


ia el eS (y—1)? sin4 bed 
ane 4, 4+ 2 (5y—3) cos? A—(y—1) cos* 0 

Sees oe ae EY 

ae 4y (1+3 cos? @)+-(y—1)? sin? 0 (5—cos? A) ss 

is (y—1)? sin4 6 , ceil? 


This last corresponds to the original Moller formula. 


§ 3. DIScUSSION 

The ratio Q(——)/Q(— <=) at the most favourable angle (9=90° in the 
centre-of-mass system) is equal to } in the limit of high energy for both e+ 
and e. As the energy of the electron decreases, this ratio decreases 
further, to zero in the low energy limit (cf. Mott and Massey 1949). This 
may be regarded as due to the effect of the Pauli principle on the prob- 
ability of the electrons being close together. The ratio for positrons tends 
to 1 at low energies. But the notion that there should be no polarization 
effect because the positron and the electron are * different ° particles would 
be quite misleading at energies of the order of mc?, where the virtual 
annihilation term is far from negligible. 

The cross sections have also been calculated by A. Bincer, according to a 
recent communication by Dr. Frauenfelder. The ratio quoted agrees with 
ours for electrons, but not in the case of positronsy. 


+ Note added in proof —These authors now obtain a slightly different formula 
(Phys. Rev., 107, 644). This is consistent with our (4) and there js no longer any 
discrepancy. 
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SUMMARY 


The deformation properties of indium antimonide have been studied 
under a variety of conditions. At room temperatures elastic behaviour 
is observed up to the fracture point, while at higher temperatures glide 
occurs on {111} planes. Impact at room temperature can cause glide 
on {110} planes. A method of observing individual dislocations is 
described, and the dislocation distributions after the various deformations 
have been studied. 


§ 1. INTRODUCTION 


THE mechanical properties of the semiconductors silicon and germanium 
have been studied by several investigators. At room temperature 
these are brittle solids showing elastic behaviour up to the fracture point 
(Pearson et al. 1957), while at higher temperatures they show plastic 
behaviour with a sharp yield point and, for stresses below this yield 
point, flow with a temperature and stress-dependent incubation period 
(Gallagher 1952). These features naturally suggest an interpretation 
in terms of a Cottrell atmosphere locking mechanism, although this is not 
the only possibility. Similar behaviour has now been observed in indium 
antimonide, a semiconductor with zinc-blende structure. In addition 
plastic deformation has been observed at room temperature under certain 
special conditions. 
§ 2. OBSERVATION OF DISLOCATIONS 


The plastic properties of crystals can usually be described in terms 
of the theory of dislocations but often more than one mechanism can be 
proposed to explain a given set of macro-properties. It is therefore useful 
if the actual numbers and positions of the dislocations can be determined 
experimentally. To this end we have devised an etch which produces on 
suitable crystal faces of InSb pits which mark the point of emergence 
of dislocations, this etch being a modification of CP4 as used by Vogel et al. 
(1953) to show dislocations in germanium. 

In this method the crystal is first ground with 850 grade Unirundum 
in aqueous Teepol suspension on a glass plate. It is then polished chem- 
ically with CP4A or with CP4A diluted with water and nitric acid, after 
which it is placed in the etchant for one minute. (The composition of 


t Communicated by the Author, 
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CP4A and of the etchant are given in an Appendix.) Best results are 
obtained if the solution is not agitated during this last process. After 
washing and drying clearly visible circular pits are observed on faces of 
the form {111}. Their visibility decreases as the angle between the face 
and {111} increases, but they may be seen as elongated ellipses on faces 
of the form {112}. It is interesting that the form {111} does not show 
these pits. This provides a rapid means of distinguishing between the 
two tetrahedral forms, a distinction necessary for the complete orientation 
of the non-centrosymmetrical class 43m but one which cannot be obtained 
by x-ray methods. (The convention adopted here as to which is the 
positive direction of the crystallographic axes conforms with that of 
Dewald (1957).) 


Fig. 1 


3x 10° 


2% 108 


Density of pits (.m-?) 


1x 10° 


0 2 A 6 8 10 12 14 
Distance across crystal diameter (mm) 


Variation of etch-pit density across the diameter of an indium antimonide 
crystal which was grown in a [111] direction. The vertical line indicates 
the statistical scatter expected. 


These etch pits are quite reproducible, in that if a surface treated in the 
above manner is ground to a depth sufficient to remove them then on re- 
etching they reappear at the same positions. By continuing the grinding 
and etching process they may be traced through a crystal and therefore they 
correspond to line imperfections. As will be shown below, their disposition 
throughout the crystal is in accord with the expected disposition of an 
array of dislocations. We therefore conclude that the position of a pit 
marks the intersection of a dislocation with the surface. The converse 
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of this proposition, namely that every dislocation site produces a pit, 
is much more difficult to prove, but experiments on bent bars have shown 
that the pit density is of the same order of magnitude as the expected 
dislocation density. 

Single crystals of InSb grown at S.E.R.L. by the Czochralski pulling 
technique with [111] as the growth axis have dislocation densities of the 
order of 2x 103 cm~?, which is of the same order as is observed in silicon 
crystals grown here by the same method, but is about three orders of 
magnitude smaller than the densities commonly observed in well-annealed 
metal single crystals. The variation of pit density across a face per- 
pendicular to the crystal growth axis for a typical crystal of InSb is shown 
in fig. 1. This is qualitatively similar to that obtained with pulled silicon 
or germanium crystals and is the result of deformation under thermal 
stresses caused by radial temperature gradients as the crystal is with- 
drawn from the melt. The small number of dislocations in this material 
helps considerably in interpreting plastic deformation experiments, 
since the density of introduced dislocations may then be three or four 
orders of magnitude greater than the density of those initially present. 


§ 3. Room TEMPERATURE FRACTURE 


Bars of InSb about 1 mm x 0-5mm x 11mm were supported at their ends 
and loaded quasi-statically at their centre. Measurements were made of 
the displacement at the centre as a function of load. Several specimens 
cleaved on (110) planes at comparatively low loads. Since cleavage 
occurred at only one place in such specimens and since the load at cleavage 
varied widely for different specimens it was assumed that this behaviour 
was due to adventitious stress concentration. Other specimens showed 
elastic behaviour up to the fracture point, at which they broke into several 
small pieces, not bounded by crystallographically important planes. 
The maximum strains at fracture were of the order of $°%, about the value 
obtained for silicon bars of similar dimensions (Pearson et al. 1957), 
Some of the fractured pieces were etched and it was found that no appre- 
ciable movement of dislocations had occurred, nor were fresh dislocations 
introduced in the bulk of the material. However, some ‘ dislocation- 
cracks ’ were seen to run from the fractured edges. This type of deforma- 
tion will be discussed below. It is evident that at room temperature 
there are considerable forces locking dislocations in fixed positions. 
This fact is of importance in those experiments in which the positions of 
dislocations are to be studied, since these positions are not changed by 
handling the crystals provided that reasonable care is taken not to subject 
them to violent stresses. 


§ 4. HicgH-TEMPERATURE FLOW 


At temperatures of 200°c or more InSb becomes plastic, and single- 
crystal bars of it can be bent through large angles, On the absolute 


Dee 


1478 J. W. Allen on the 


scale 200°c is about two-thirds of the melting temperature, and approxi- 
mately the same ratio of temperature required for plastic flow to melting 
temperature is observed in silicon and germanium. In atypical experiment 
a bar 0-5 mmx 1 mmx 15 mm is supported at its ends and a fixed load is 
applied at its centre, this load being of the order of 50 gm wt, i.e. about 
a tenth of that which would cause fracture at room temperature. This 
system is placed in an oven containing a window and measurements of 
displacement at the centre as a function of time are taken. The variation 
obtained is shown diagrammatically in fig. 2; it is very similar to that 
observed in germanium (Gallagher 1952). This graph has three different 
regions corresponding to an incubation period, plastic flow, and work 


hardening. 
Fig. 2 


cn a 


Strain ———> 


Time ———> 


Strain—time curve for an InSb bar bent at an elevated temperature. 


During the incubation period AB there is no appreciable deformation. 
The length of this period decreases as the stress or the temperature 
increases, and for the system described above is of the order of twenty 
minutes. It does not appear to be changed by preheating, i.e. by holding 
the crystal at the temperature of the experiment for some time without 
aload. Seitz (1952) has put forward an interpretation of the incubation 
in terms of the diffusion of locking impurities away from dislocation 
sources, while Pearson et al. (1957) suggest the Peierls force as being the 
source of locking. Experiments are in hand to determine, by the etch-pit 
technique, the number and disposition of dislocations introduced under 
various conditions of deformation and in this manner to unravel the 
processes involved. 

The period of incubation is followed by one of plastic flow, BC. After 
deformation slip lines can be seen on the crystal which serve to establish 
the glide planes as being {111}, the planes of closest packing. The 
visibility of the slip lines on different faces of the bar depends on the 
crystallographic orientation and on their orientation with respect to the 
bend axis. ‘This effect may be used to determine the slip direction, and 
this appears to be [110], the direction of the shortest translation vector 
in the zince-blende lattice. On etching, the dislocations introduced may be 
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seen. Their general disposition is shown in fig. 3, Pl. 60, which is a (IIT) 
face of InSb viewed along the bend axis. It will be seen that dislocations 
approach the neutral plane from both sides ; they are generated mainly 
by sources at or near the upper and lower surfaces which are, of course, 
regions of maximum stress. While many of the dislocations are con- 
centrated into relatively few slip lines an appreciable number are not. 
Another point of interest is that the arrays are not strictly linear, the 
deviations from linearity being of the order of a few microns. This is not 
due to interaction with the surface during slip, since by grinding off the 
surface the same effect may be seen to have taken place within the bulk 
of the crystal. Nor is climb out of the glide pane by thermal diffusion 
a likely explanation since this specimen was bent under conditions such 
that the incubation period was very short, and was air-quenched 
immediately after deformation so that the total time from applying a load 
to the crystal to its returning to room temperature was about three 
minutes. At 260°c, the temperature used in this particular experiment, 
the distance through which a dislocation will diffuse thermally in three 
minutes is certainly less than the distances with which we are here con- 
cerned. One possible explanation is that there is cooperative slip on 
parallel glide planes, although the scale here is much larger than is observed 
in multiple slip in aluminium (Heidenreich and Shockley 1948). Some 
credence is lent to this hypothesis by fig. 4, Pl. 60 which shows part of 
fig. 3, Pl. 60 at greater magnification, and in which some sort of cooperative 
slip may be seen. Another possibility is that the obstacles to movement 
presented by impurities are sufficient to make short lengths of dislocation 
climb out of their glide plane into another parallel one, although the 
magnitudes of the deviations are rather large to be accounted for by such 
a process. It should also be noted that the dislocations do not lie in the 
equilibrium positions calculated by Eshelby eé¢ a/. (1951), which implies 
that there is considerable resistance to their motion even in the glide 
plane. 

Plastic flow continues for some time, after which work-hardening sets 
in, as indicated by CD in fig. 2. If neither temperature nor load are too 
great this part of the strain-time curve becomes essentially flat. If now 
either temperature or load is increased further deformation takes place, 
after which the curve flattens off at a higher value of the strain. The 
amount of deformation before work-hardening also depends on orienta- 
tion : if the orientation is chosen so that slip on one set of glide planes is 
preferred the section BC of fig. 2 is considerably extended. By etching 
to reveal dislocations it can be shown that work-hardening is associated 
with slip on two or more intersecting glide planes. An example of such 
an intersection is illustrated in fig. 5, Pl. 60; a short etching time has 
been used in order to show more clearly details of the dislocation positions 
at the crossing point. Some ‘ piling-up’ can be seen. This not so pro- 
nounced as in a metal, for in InSb under these conditions the resistance 
to glide is high in any case so the intersecting glide plane does not increase 
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this resistance by a large factor. Observations of these intersections | 
also provide evidence that the dislocations themselves must to some 
extent have crossed during slip, and the jogs thereby produced presumably 
contribute to work-hardening. 


§ 5. RooM-TEMPERATURE FLOW 


Tt has been shown above that quasi-static bending of InSb at room 
temperature causes fracture before an elastic limit is reached. However, 
if InSb crystals are subjected to sudden local impact ‘ dislocation-cracks * 
may be formed. These are observed as line markings on the surface and 
have the appearance of fine cracks under low magnification, while under 
high magnification after etching individual dislocations may be resolved 
along their length, or part of it. In general these ‘ dislocation-cracks * 
are not straight, although they show a strong preference to lie in {110} 
planes. The dislocation density is not uniform along them, although 
short sections may show uniform spacing. They rarely run right across a 
crystal, and they are sometimes the continuation of an actual crack. 
Normally two or three more widely spaced dislocations may be seen along 
the continuation of the line some microns beyond its termination. Pre- 
liminary experiments suggest that the most favourable conditions for 
their formation are that a stress of the order of the fracture stress should 
be applied suddenly to a small region of the crystal, the rest of the crystal 
being comparatively stress-free. The dislocations are then propagated 
from some source and run into the crystal until the applied stress is 
removed or until they reach a region where the internal stress is in- 
sufficient to keep them in motion. The fact that on removal of the 
applied stress the dislocations do not run out again, and also that their 
spacing in the slip plane is very non-uniform, shows that there is a strong 
locking force which in this case cannot be a Cottrell atmosphere, since 


the times involved are far too short for the formation of such an 
atmosphere. 


§ 6. CONCLUSIONS 

Indium antimonide is a convenient substance for studying the pro- 
duction and properties of arrays of dislocations and for comparing these 
with theoretical ideas ; experimentally, dislocations may be introduced 
under controlled conditions, they retain their positions and these positions 
are easily determined, and theoretically, the structure is not of intractable 
complexity. The low density of dislocations in our grown crystals is 
also of considerable advantage. As an example, a sample of reasonable 
size may only contain of the order of ten dislocations, and so nucleation 
of further dislocations at the surface is not masked by nucleation at 
Frank—Read or other internal sources. The mechanical properties of 
InSb are similar in many respects to those of a face-centred cubic metal 
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and to those of silicon and germanium, although it shows some character- 
istic differences; we may therefore hope that investigations of these 
properties in InSb will not only be of intrinsic interest but will shed light 
on the behaviour of these other materials. 
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APPENDIX 
CP4A polishing etch 111 etch to show dislocations 
5 parts HNO, aq. 1 part CP4A 
3 parts HF aq. 1 part H,O dist. 


3 parts CH,COOH glacial 1 part CH,COOH glacial 
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Electron Microscopic Observation of Periodic Structures 
below 104. 


By G. A. Basserr and J. W. MENTER 
Tube Investments Research Laboratories, Hinxton Hall, Cambridge 


[Received August 19, 1957] 


In his initial studies on the direct resolution of lattice planes in crystals 
by electron microscopy one of the authors (Menter 1956) worked with 
metal phthalocyanines, principally because of their relatively large 
lattice parameters (e.g. d(g3) in platinum phthalocyanine=11-94 A). 
The remarkable clarity of the images obtained suggested that it should be 
possible to resolve more closely spaced lattice planes. This was supported 
by an elementary consideration of the effect of the lens aberrations on 
the resolution of periodic structures which showed that these are less 
deleterious than with non-periodic structures. 

However, for the practical realization of improved resolution a number 
of exacting conditions must be fulfilled by the specimen. It must be 
reasonably stable in vacuum and able to withstand electron bombard- 
ment without decomposition. The habit of the crystal must be such as 
to present the lattice planes it is hoped to resolve closely parallel to the 
electron beam when the crystal is resting on a support film. The crystal 
must be thin (a few hundred A) in the direction of the electron beam. 
The number of crystals fulfiling these conditions with planar spacings 
between 5 and 10 A is not great. At the present time the most successful 
results have been obtained with crystals of molybdenum trioxide in which 
the (020) planes (d=6-93 A) have been clearly resolved. Molybdenum 
trioxide has perfect cleavages on (010) and (100) and a less perfect but 
distinct cleavage on (001) so that there is a reasonable probability of 
obtaining suitably oriented fragments. Crystals have been prepared by 
two methods. (1) A large crystal kindly supplied to us by Mr. H. P. 
Rooksby of G.E.C. Research Laboratories, Wembley was broken down 
mechanically by grinding. (2) Small crystals were prepared by thermal 
decomposition of ammonium molybdate followed by resublimation. In 
both cases the final preparation was suspended in alcohol for deposition 
on electron stabilized nitrocellulose support films containing many holes. 

A considerable number of micrographs have been taken using a Siemens 
Elmiskop I with an instrumental magnification of x 160 000, accelerating 
voltage 80 kv, objective aperture 50 p. They were recorded on Ilford 
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Contrasty Special Lantern plates developed in [ford developer ID36. 
About 10% of the plates show fragments of crystal in which the lattice 
has been resolved. The spacing measured from the micrographs lies 
within --10% of the x-ray value quoted above. Due to the decomposi- 
tion of the crystals in the microscope it has not been possible to record 
an electron diffraction pattern and a micrograph showing the resolved 
lattice from the same crystal in order to establish the orientation 
unambiguously. Assuming the habit surfaces are as defined by the 
cleavages given above then (020) planes may be resolved in crystals lying 
either on (001) or (100). 

The clearest micrograph obtained in this investigation is shown in 
fig. 1, Pl. 61. The crystal is attached to the edge of a larger crystal itself 
spanning a hole in the support film. The background structure of the 
latter has been avoided in this way but the quality of the image is still 
seriously impaired by the structure of the contaminant carbon film C 
which envelops the crystal rapidly in the course of examination. A pro- 
jection of the crystal structure on (001) as determined by Wooster (1931) 
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Schematic projection of molybdenum trioxide structure on (001), showing 
positions of molybdenum atoms only (after Wooster 1931). Same 
scale as fig. 1, Pl. 61. 


showing only the positions of the molybdenum atoms is shown in fig. 2 
to the same scale as the electron micrograph. It is clear from this how 
far the electron microscope image still is from representing the true 
structure of the crystal due to the limited number of diffraction spectra 
contributing to the image. 

It is considered that this resolution of a singly periodic structure of 
6-9 A by no means represents the limit of what is possible with existing 
microscopes. In the course of other experiments in our laboratory, 
moiré patterns formed by (422) planes in a gold/nickel specimen have 
been clearly resolved with a spacing of 5-3 A} (Pashley et al. 1957). 
Furthermore doubly periodic moiré patterns with a spacing of 9-2 At 
from the (220) planes of the same combination of metals have also been 
resolved showing that lens astigmatism is not yet a limiting factor. The 
clarity of these micrographs suggests that doubly periodic structures 
below 9 A should be readily resolvable. Since the mechanism of forma- 
tion of the moiré image by recombination of diffracted spectra is essen- 
tially similar to that of a directly resolved lattice, we believe that with 
suitable crystals direct observation is possible at least down to 5 A. 


+ These micrographs are in the course of publication elsewhere. 
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Absorption lines of MnCl,.4H,O at 78°K. 
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(a) (d) (c) 


Laue transmission photographs, quartz single crystals oriented with incident 
X-ray beam parallel to c-axis. (a) Before irradiation, (b) after 
7X 101% neutrons/em?, (c) after 1-2 x 102° neutrons/cm2. 
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(a) (0) (c) 


Precession photographs, zero level (hkO) projection, quartz single crystal. 
(a) Before irradiation, (6) after 7x 10!% neutrons/cm?, (c) after 
1-2 « 107° neutrons/em?,. 


Fig. 10 


(a) (b) 


X-ray transmission scattering patterns from neutron irradiated quartz single 
erystals. Monochromatic Mo,;,radiation. (a) After 7 x 10" neutrons/cm’, 
(6) after 1-2 x 107° neutrons/cm?. 
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Dislocations in a bar of InSb bent through a small angle. (111) face. 47. 


An enlargement of part of fig. 3 showing slip on neighbouring glide 
planes. x93. 


Fig. 5 


Dislocations in intersecting glide planes. (111) face of InSb. 180 
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Fig. 1 
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Electron micrograph of molybdenum trioxide crystal showing (020) planes 
resolved. (d,g09) 6:93 A.) Magnification x3 000 000. 
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